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SYNOPSIS 

A STUDY OF OVERVOLTAGE PROBLEMS IN EHV SYSTEMS 


by 

N.S. Saxena 
Ph.D. 

Department of Electrical Engineering 
Indian Institute of Technology, Kanpur 
July, 1981 


With the increase in power demand, the present trend 
is to transmit large amount of power over long distance 
lines at extra high voltage (EHV). In EHV systems, the 
insulation level is decided by overvoltages caused due to 
switching and/or faults rather than lightening. A poor choice 
of insulation level will either increase the cost of the 
system or will increase risk of failures due to faults 
initiated by flashovers. Therefore a proper insulation level 
should be chosen which requires the knowledge of the nature 
and magnitudes of various overvoltages caused by different 
system conditions. Since the magnitude and time duration of 
various overvoltages is different, therefore, they stress the 
insulation in different manner. 

There are three distinct periods during which the 
overvoltages occur, viz., 
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a) A surge period during which the travelling wave effects 
predominate and in which the line is represented by its 
surge characteristics. The transient (surge) overvoltage 
are caused by line energisation/re-energisation or fault 
initiation. 

b) A dynamic period which is transitional between the surge 
period and the steady state and is characterized by the 
voltage variation contained in an envelope which varies 
aperiodically with time. 

c) A steady state period during which the voltage is 
periodic but usually is distorted due to various harmo- 
nic components. When the power system is terminated by 

a transformer, steady state overvoltage may occur at the 
transformer terminals due to switching operation. The 
transformer non-linear characteristic plays a dominant 
role in this phenomena and the system is said to be in 
a ’ ferro-resonant * mode. 

The objective of this thesis is to study all the 
above mentioned overvoltages. Suitable techniques for digi- 
tal simulation have been evolved for their evaluation which 
consider more detailed system representation than reported 
before in the literature. 

Steady state overvoltages are studied first. The 
field investigation of an electromagnetic voltage transfor- 
mer (EMVT) failure while being de-energized by a circuit 
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breaker (CB) fitted with grading capacitors across its 
contacts is described. This phenomena has been analyzed 
using the principle of ’Harmonic Balance' and it is shown 
that the failure was due to f erro-resonance phenomena. Also 
the analysis for the pre-steady state period is carried out 
using state space techniques. This is to give insight into 
initial conditions which might give rise to the steady state 
overvoltages. Based upon the data collected for several 
makes of EMVTs and CBs fitted with grading capacitors, reme- 
dial measures have been suggested to avoid such incidents. 

Dynamic overvoltages following load rejection have 
been considered next. The methods reported in the literature 
for the evaluation of such overvoltages are approximate. 
Usually the digital simulation is employed to predict the 
variation in magnitude and frequency of the internal voltage 
of the machine and this information is used subsequently 
on the transient network analyzer (TNA) to predict ferro- 
resonant overvoltages. The digital simulation ignores the 
network transients but can predict the phenomena of self 
excitation. In this thesis, this two stage computation is 
combined in a single digital simulation program, which 
incorporates not only the network transients on a 3-phase 
basis but also the synchronous machine representation in 
detail. The prime mover-governor is also considered along 
with the exciter and voltage regulator. The nonlinearities 
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due to transformer and machine saturation are included. It 
is belived that this is a more accurate approach, and can be 
used to determine if the phenomena of self excitation and/or 
ferror-resonance can occur in a single simulation. This is 
illustrated with the help of the case studies of both hydro 
and thermal machines feeding a load through a long EHV line. 
The effectiveness of the voltage regulator in controlling 
the overvoltages due to f erro-resonance is also demonstrated. 

Transient overvoltages are considered next. The 
method of characteristics as given by Dommel is recognized 
to be most efficient for the calculation of the transmission 
line transients, since the equations involved are algebraic 
rather than differential. However this method requires 
iterative solution at each time step if nonlinearities are 
present in the network. The state space techniques are more 
efficient in handling nonlinearities. In this thesis the 
alternative of (l) employment of state space techniques for 
both the network and the transmission line and (2) the 
combined use of state space techniques for the network and 
the method of characteristics for the transmission lines, 
have been examined. It is found out that the second approach 
is usually superior in terms of computation time. 

The transient overvoltages caused by fault initiation 
are also considered. The transmission line model using the 
method of characteristics due to Dommel is extended to 
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simulate the initiation of a fault at any point of the line. 
The generator dynamics is represented, and the results are 
compared when the generator is represented by a constant 
voltage source. It is demonstrated that the overvoltages are 
greater if the generator dynamics is included in the simula- 
tion. It is found that the opening of the CB at the receiving 
end in trying to clear the fault aggravates the overvoltage 
problem and it can be controlled through the use of CB 
opening resistors. It is also found in this case that the 
overvoltages are higher at receiving end when the fault is 
closer to the receiving end. 

Finally, the thesis concludes with the highlights of 
the work reported in the area of overvoltages in EHV systems 
and projects the scope for further work. 



CHAPTER - I 


INTRODUCTION 

1.1. NATURE OF OVERVOLTAGES IN EHV SYSTEMS 

With the increase in power demand and the economic 
necessities to locate generating stations in remote areas, 
the present trend is to transmit the bulk power over long 
distance extra high voltage (EHV) lines. In an EHV system, 
the insulation level is determined by overvoltages caused due 
to switching and/or faults rather than lightening. A poor 
choice of insulation level will either increase the cost of 
the system or will increase the risk of ^failures due to 
faults initiated by flashovers. Therefore a proper insulation 
level should be chosen, which requires the knowledge of the 
nature and magnitude of various overvoltages caused by diff- 
erent system conditions. 

Following a disturbance, there are three distinct 
periods during which the overvoltages occur, viz., 

a) A surge period during which the travelling wave 

effects predominate and in which the line is represented by 
its surge characteristics ; The transient (surge) overvoltages 
are caused by sudden switching or fault inception. 
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b) A dynamic period, which is transitional between the 
transient a nd the steady state and is characterise^ by the 
voltage variation contained in an envelope which varies 
aperiodically with time. 

c) A steady state period during which the voltage is 
periodic but usually is distorted due to various harmonic 
components* When a power system is terminated by a transformer, 
steady state overvoltages may occur at the transformer ter- 
minals due to switching operation. The transformer .non-linear 
characteristic plays a dominant role in this phenomena and 
the system is said to be in a ' ferroresonant mode'. 

We give below a review of the significant work 
already done in the above areas along with the modifications 
and additions implemented in this thesis. The steady state 
overvoltages have been studied first, followed by dynamic 
and transient overvoltages respectively. 

1.2. STEADY STATE OVERVOLTAGES 

As already mentioned in Sec. 1.1, steady state 
overvoltages may occur when a system is terminated by a 
transformer. The combination of the circuit parameters may 
be such that upon switching, the phenomena of ferroresonance 
takes place. This is particularly true when air circuit 
breakers with grading capacitors are used for switching 
operation. Such overvoltages have been reported to result 
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in damage both to the electromagnetic voltage transformer 
(EMVT) and the power transformer [24,25,32,33,41]. 

The problem ha$ been analysed in the pa$t either 
by using (a) state space technique [l8,42]or (b) dual input 
describing function technique [13,35], 

The non-linearity arising due to the transformer 
saturation is approximated either by an exponential fun- 
ction [18] or by a piece-wise linear model [42]. 

In state space technique, a set of first order non- 
linear differential equations are solved on the digital 
computer. If the solution is carried out for a sufficiently 
long time, the steady state solution can be obtained. The 
steady state solution consists of the fundamental and higher 
harmonics. However, this method does not provide any direct 
information regarding the harmonic components. The nature of 
the solution also depends on the initial conditions as it 
may not always lead to overvoltages. If one wants to examine 
the effect of several system parameters, then the simulation 
of all possible initial conditions with different combination 
of system parameters becomes a laborious if not an impossible 
task. 

Swift [13] has employed dual input describing 
function technique which is popular in control systems 
theory. In this technique a small perturbation is considered 
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and the describing function is obtained assuming the non- 
linearity to be a third degree polynomial. Only the funda- 
mental frequency component of the solution is obtained since 
the overall system is assumed to be a low pass filter. The 
condition of the instability i.e. ferroresonance jump is 
obtained by calculating the critical magnitude of disturbance 
for which the Nyquist diagram encloses the minus one point. 

Prusty and Sanyal [35] generalised the Swift's technique 

is h, 

so that the non-linearity can be approximated by an n 
degree polynomial. The drawback of the Swift’s technique is 
that the magnitude of the overvoltages predicted does not 
include the effect of harmonics. Also it gives the critical 
value of the disturbance for which the overvoltages will 
occur* but no information is available about the system ini- 
tial conditions which lead to the f erroresonant mode. 

In this thesis, the principle of 'Harmonic Balance' 
[2,3] if used to determine the steady state overvoltages. 

The method gives the magnitudes of both the fundamental and 
harmonic components separately. It gives various possible 
modes of the operation and their stability can also be exa- 
mined. The proposed method is shown to give excellent agree- 
ment between the results obtained fyom the analysis and 
actual field test of a practical system* By collecting the 
data for several makes of circuit breaker and electromagnetic 
voltage trans forme r* , we have carried out detailed investigation 
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in an attempt to suggest remedial measures to prevent the 
occurrence of such overvoltages. It is shown that if the 
transformer saturation characteristic can be controlled to 
a reasonable extent, such overvoltages do not occur. 

1.3. DYNAMIC OVERVOLTAGES 

Sometimes situations may arise that -due to sudden load 
rejection, an unloaded EHV line is left connected to over- 
speeding generator(s) . This situation can lead to excessive 
overvoltages in the system. Essentially such overvoltages 
are of fundamental frequency with harmonic components and 
may persist for a few seconds. They can become dangerous 
if the generator self-excitation takes place due to excessive 
line charging and/or when ferroresonance takes place due to 
the transformer saturation. The magnitude of such overvoltages 
is influenced by the generator dynamics, excitation system, 
prime mover governor, line and transformer constants and 
saturation characteristics of the transformer, and generator. 

I 

Therefore, all the components described above need to be 
modelled accurately for the correct evaluation of dynamic 
overvoltages . 

Demello -et -al [4] employed an electronic differential 
analyzer for the determination of dynamic overvoltages. They 
represented the generator in details including amortisseur 
circuits, using the Park ’ s equations. The excitation system 
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was represented in detail while the prime mover governor 
was represented approximately by a speed versus time curve. 
Besides the above approximation, the electrical network 
was represented on a steady state basis. While the method 
could predict the generator self-excitation, it was not 
capable of predicting ferroresonance, since the network 
transients were ignored. 

Aggarwal et al [26,27] employed state space approach 
to study this problem. The system was represented on a single 
phase basis, which is not adequate when the non-linear elements 
are considered. The rise in frequency and generator internal 
voltage was represented by an approximate method. The non- 
linearity due to transformer saturation was considered by 
piecewise linear approximation. While the method was capable 
of predicting the ferroresonance, it could not predict the 
generator self-excitation. 

Clerici and Didriksen [23] employed digital simula- 
tion to predict variation of the magnitude and frequency of 
the generator internal voltages and this information was 
subsequently used on a 'Transient Network Analyzer' (TNA) 
to determine dynamic overvoltages. The digital simulation, 
where the network transients were ignored, could predict 
the generator self-excitation and the TNA could simulate the 
network transients including ferroresonance. However, the 
process required a two stage computation. 
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In this thesis we have employed digital simulation 
technique, which incorporates not only the network transients 
on three phase basis but also the generator in complete 
details including amortisseur circuits. Its excitation system 
and prime mover governor are also represented in detail* 

The non-linearities due to transformer and generator saturation 
are also included. The synchronous generator is modelled by 
a dependent current source in parallel with a constant indu- 
ctance as given by Ramshaw and Padiyar [28]. The time 
varying coupling in the machine inductance matrix is repla- 
ced by a dependent current source, The concept of dummy rotor I 

I 

coils is used to eliminate dynamic saliency and thereby 
avoiding the inversion of time * varying matrix at each step 
of integration. 

It is believed that this is a more accurate approach 
for the determination of dynamic overvoltages and can predict 
the phenomena of self-excitation and/or ferroresonance in a 
single stage simulation. This is illustrated with case studies 
for both hydro and thermal generators feeding a load through 
a long EHV line. The effect of the excitation system in 
controlling the overvoltages due to ferroresonance and self- 
excitation is also demonstrated. 
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1.4. TRANSIENT OVERVOLTAGES - SYSTEM MODELLING 

The transient overvoltages can arise in an EHV system 
either due to internal causes such as switching or fault 
inception or due to external causes such as lightening die- 
charges. It is a well known fact that the EHV systems are 
designed from the consideration of overvoltages due to 
switching [5, 14] or fault inception [15,16,17]. The calcula- 
tion of these overvoltages requires the representation of 
the transmission line on a distributed parameter basis and 
hence is quite complex. There are several method [6, 7, 8, 9, 

10,, 11, 12] for the transmission line modelling, among them 
the method of characteristics suggested by Dommel [11,21,31] 
is perhaps most widely used. Using this model for the trans- 
mission line, the resulting equivalent circuit for the EHV 
system is purely resistive and hence can be easily solved. 
Dommel [2l] proposed 'Compensation Method' when non-linearities 
due to the transformer saturation etc. have to be considered. 
However, the method becomes very slow due to its iterative 
nature when there are several non-linearities. 

The non-linearities can be better handled with the 
use of state space techniques. Therefore, in this thesis, 
two alternatives of the system modelling have been considered 
viz., (a) to interface the transmission line model based on 
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the method of characteristics with the state space model for 
the rest of the system and (b) to employ the state space 
model for the transmission line as well along with the rest 
of the system. It is shown that the Alternative 'a' is superior 
to Alternative *b' in terms of the computation time. 

1.5. transient overvoltages due to faults 

The transient overvoltages which are considered impor- 
tant for the insulation design for EHV systems are due to line 
energisation/re-energisation [5,14,27]. With further increase 
in the system voltage, the safety margins are reduced due to 
economic reasons. This requires that the line energisation/ 
re-energisation voltage should be controlled to a lower value 
which is achieved mainly by circuit breaker closing resistors. 
In such cases the overvoltages due to fault inception and 
subsequent load throw-off may become more important [16,23]. 

The study of such overvoltages requires detailed generator 
representation especially for the salient pole machine [lj. 

In earlier TNA studies [16,22,23], the detailed generator 
representation was not possible. A method to interface the 
detailed generator model with the electromagnetic transients 
program was given by Brandwajn and Dommel [39]. The method 
employed Park's equations for the simulation of the generator 
and a few approximations were made to avoid the re- 
tringularisation of the network matrix at each integration 
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step due to the time dependent generator inductance matrix 
introduced by dynamic saliency. 

The generator representation employed in this thesis 
has been discussed in Sec. 1.3. It may be noted that this 
representation gives constant inductance matrix without 
making any approximation. Based upon the model derived from 
considerations discussed in Secs. 1.4, the study of over- 
voltages due to fault inception and subsequent load throw-off 
have been carried out. The effect of the detailed generator 
representation, number of generators and the location of the 
fault have been studied. 

1.6. OUTLINE OF THE THESIS 

Chapter-II deals with the determination of steady 
state overvoltages. The failure of an electromagnetic voltage 
transformer in a practical system while being de-energised by 
a circuit breaker fitted with grading capacitors has been 
studied. The field investigations subsequent to the damage 
and the analytical study performed using the principle of 
'Harmonic Balance' have been presented. The field and ana- 
lytical results obtained are shown to be in close agreement. 
The analysis of the pre-steady state period is also carried 
out to get an insight into the initial conditions which give 
rise to the steady state overvoltages. This analysis is done 
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using the state space techniques. Based upon the data 
collected for several makes of electromagnetic voltage trans- 
formers and circuit breakers fitted with grading capacitors’, 
remedial measures have been suggested to avoid such incidents. 

Chapter - III deals with the determination of dynamic 
overvoltages following load rejection. A digital simulation 
technique for the evaluation of such overvoltages has been 
presented, which not only accounts for the network transients 
on three phase basis but also includes the detailed generator 
representation along with its excitation system and prime 
mover governor. The non-linearities due to transformer and 
generator saturation are also considered. The simulation 
technique is capable of predicting the self-excitation and/or 
f erroresonance in a single simulation. Both the hydro and 
thermal systems are considered and the effectiveness of an 
excitation system in controlling the overvoltages due to self 
excitation and/or ferroresonance is illustrated. 

Chapter IV deals with some modelling aspects for 
the determination of transient overvoltages. The method of 
characteristics is more popular in determination of transient 
overvoltages. However the EHV system has several non-linear 

elements, the state space technique may be more convenient. 

Thus, two alternatives have been considered viz., (a) to 
interface the method of characteristics for the transmission 
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line representation with the state space model for the rest 
of the EHV system and (b) to use the state space model for 
the entire EHV system including the transmission line. Both 
the alternatives have been studied and it is concluded that 
Alternative (a) is superior to Alternative (b) . 

Chapter V deals with the transient overvoltages 
caused by fault inception and also subsequent load rejection. 
The method of characteristics is extended to simulate the 
fault at any point on the transmission line. The fault 
inception overvoltage are studied for both detailed and 
simplified generator representation. The effect of the number 
of generators operating in parallel and fault location is 
also studied. Subsequent to the above studies, the over- 
voltages on load throw-off» while the fault (at various 
location on the line) is still present, are obtained. The 
effectiveness of the circuit breaker opening resistors in 
controlling such overvoltages is demonstrated. 

Chapter VI concludes with a review of the significant 
work done in this thesis and scope for further work. 
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CHAPTER II 

STEADY STATE OVERVOLTAGES 

2.1. INTRODUCTION 

Steady state overvoltages may occur jWhen a system 
is terminated by a transformer. The combination of the 
circuit parameters may be such that upon switching the 
phenomena of f erro-resonance takes place. This is particu- 
larly true when air circuit breakers with grading capacitors 
are used for switching operation. Such overvoltages have 
been reported to result in damage both to the electromag- 
netic voltage transformer (EMVT) and power transformer 
[24,25] and, thus, have been the cause of concern. 

The main approaches adopted so far in the literature 
to study this problem have already been reviewed in 
Chapter I, where the shortcomings of the state space and 
describing function techniques for steady state solutions 
have been given. 

In this thesis, we have adopted the principle of 
'Harmonic Balance' [2,3] for studying steady state over- 
voltages. We were fortunate enough to analyse the failure 
of a 220 KV electromagnetic voltage transformer (EMVT) which 
occurred in the U.P. State Electricity Board (UPSEB) system 
in India. Subsequent to the damage, field investigations 
were carried out to simulate what actually happened and is 
given in Section 2.2. 
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2.2. FIELD INVESTIGATIONS 

A 220 KV bus EMVT installed at a 220 KV substation 
(whose single line diagram is given in Fig. 2.1) got damaged 
when the 220 KV bus was to be completely de-energised from 
all sources of supply by opening C, the last circuit breaker 
(CB) (associated bus isolators remaining closed) when other 
CBs (A and B) had already been opened. 

The damaged EMVT was replaced by an EMVT of the 
same make and design and the sequences which resulted in 
damage to the earlier EMVT were rj-enacted and it was obser- 
ved that when the last CB was opened to disconnect the 
source of supply, the voltage instead of falling to an 
anticipated low value, some times shot to a high value of 
380 KV. 

The observed voltage waveforms before and after 
opening of the CB are shown in Figs. 2.2A, 2.2B and 2.2C. 
From Fig. 2.2A, it would be seen that on opening the CB 
at point A, there is a transient period before the voltage 
settled to a steady state value of 380 KV RMS from point B 
onwards* The transient period was found to vary as seen 
from Figs. 2.2A and 2.2B. Figure 2.2C shows a case when 
the voltage settled to a low value of 70 KV RMS. These 
results were also repeated with an additional loading of 
9 Ohm resistor connected across EMVT secondary to dampen 
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220 KV LINE I 220KV LINE H 



1 . X DEPICTS AIR BLAST CIRCUIT BREAKER 

2. DEPICTS ISOLATORS 

3. } DEPICTS EMVT (ONE ON EACH PHASE) 

4. ALL ISOLATORS WERE IN CLOSED POSITION 

FIG. 2.1 SIMPLIFIED SINGLE LINE DIAGRAM OF 
220 KV SUB-STATION- 
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<C> 

2.2 VOLTAGE WAVE FORMS BEFORE AND AFTER 
CB OP E NING . 
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the oscillations but no significant change in results was 
observed. 

The v-i characteristic of the EMVT as shown in 
Fig. 2.3A was obtained from the field test. The piece-wise 
linear approximation was used later in the analytical study. 
It may be noted that the slope is not known from the 
field data and, therefore ,is a variable. . 

The oscillograms of the exciting current upto the 
operating voltage range when excited from the secondary and 
the primary side are given in Figs. 2.4A and 2.4B respect- 
ively. Sine v-i characteristic of the ' EMVT is nonlinear, 
these waveforms are rich in harmonic content. 

2.3. ANALYTICAL STUDY 

The equivalent circuit representation of the power 
system as seen from the 220 KV bus is given in Fig. 2.5. 

The various quantities shown in Fig. 2.5 are as 
follows : 

= Source voltage 

R and X = Source resistance and reactance respectively 

y y 

C g = Circuit breaker grading capacitance 

R = Resistance to account for losses in the EMVT 
and connected apparatus at the bus plus the 
bus loss itself. 

L = EMVT nonlinear inductance. 
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FK3.2.3A v-i CHARACTERISTIC OF EMVT. 



FIG.2.3B PECEWISE LINEAR APPROXIMATION OF 
THE EMVT v-i CHARACTERISTIC • 



FIG. 2.4 A EXCITATION CURRENT WAVE SHAPE 
WITH SECONDARY ENERGISED. 



FIG. 2. 40 EXCITATION CURRENT WAVE SHAPE WITH 
PRIMARY ENERGISED. 
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c g Cg 



FIG.2.5 EQUIVALENT CIRCUIT REPRESENTATION 
OF THE POWER SYSTEM AS SEEN 
FROM 220 KV BUS. 



FIG. 2. 6 EQUIVALENT ORCUIT OF THE SYSTEM 
NEGLECTING SOURCE IMPEDANCE - 
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If we neglect the source impedance, the equivalent 
circuit of the system is as shown in Fig. 2. 6. The system to 
the left of the terminals 'a' and ’b' is represented by its 
Thevenin's equivalent and the resultant system is as shown 
in Fig. 2.7. 


The objective is to find the voltage across L of the 
EMVT when the CB is opened. This can be done by solving for 
flux 'cp' developed in the EMVT core since the voltage is the 
derivative of flux. Therefore, the v-i characteristic of 
the EMVT given in Fig. 2.3B is converted to a <p -i chara- 
cteristic (Fig. 2.8) as seen from the primary side. 


12 + - 
dt C 


th 


f i r dt 


'th 


( 2 . 1 ) 


i R + f ( <p) 


( 2 . 2 ) 


•R 


R. i 


R 


(2.3) 


If we take 


'th 


Eth sin ( Art + &) 


(2.4) 


Then equation (2.1) can be written after differentiation as 

* 2 , 

L a. Tr 

dt 


+ k f(<p) = G cos (u+6 ) 


(2.5) 


dt 


'th 


where , 



k 


( 2 . 6 ) 



FLUX <MWb) 
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PR1MARRY CURRENT Ip(mA) 


FIG. 2.8 <M CHARACTERISTIC OF VT . 
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G = a) e.^ (2,7) 

u = cut ( 2 . 8 ) 

Both steady state and pre-steady state period ana- 
lysis of the system were done. The steady state analysis is 
given in Section 2.3.1 and the pre-steady state period 
analysis is given in Section 2.4 

2.3.1. Steady state analysis 

Vie wish to solve eqn.(2.5) for the steady state 
flux , as the derivative of flux <p is the desired volt- 
age across EMVT. Since the equation is nonlinear, the flux <p 
and the magnetising current f(<p) can both have harmonic 
components besides the fundamental. The principle of 
•harmonic balance' [3j is used to solve for the flux 9 
For the time being, we will concentrate on the fundamental 
component and determine the harmonic components later, we 
can assume 9 of the form 

9 = X cos u (2.9) 

As seen from eqns. (2.5) and (2.9), 6 is the phase 

difference between the forcing function g and the flux 9 . 

Substituting eqn. (2.9) in eqn. (2.5), we get 
X cos u - k*>.X.sin u + f (<p ) / = G c cos u + G g sin u 


( 2 . 10 ) 
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where, 

G = G cos 5 (2.11) 

c u 

G = G sin 6 (2.12) 

O 

The term f(X cos u) in eqn.(2.10) is a periodic 
function of time and hence it can be resolved by Fourier 
series and written as 

f (9 ) = A + 2 [A sin(nu) + B cos(nu)] (2.13) 

n n n 

Since flux 9 is an even function of time, therefore, 
f ( 9 ) will also be an even function of the time. This can be 
seen from Fig. A-l in the Appendix. Therefore, the sine 
components of f( 9 ) will be zero. Further, since the steady 
state f(<p) does not have dc component, A q is also zero. As 
we are for the moment considering only the fundamental comp- 
onent, f(<p) can be written as 

f( 9 ) = B^ cos u (2.14) 

An expression for has been derived in Appendix A. 
As seen from eqn.(A. 8 ) and Fig. A-l it is a function of X 
(peak amplitude of flux 9 ), break point fluxes 9 -^ and <p 2 
and slopes L^, and L 3 (see Fig. 2.8). Substituting eqn. 
(2.14) in eqn. (2.10) and collecting the coefficients of 
cos u and sin u separately and equating them we get, 

-^X + B x /C th » G<, (2.15) 


-ku X 


(2.16) 
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where, 

G 2 = G 2 + G 2 (2.17) 

c s 

Equations (2.15) and (2.17) can be solved fox flux 
amplitude X iteratively by 'Secant* method (Appendix B). 

The steps to be followed are as follows t 

1. Assume two values for X viz. X Q and X^ . 

2. Obtain the corresponding values of B, from 
eqn, (A*8). 

3. Obtain the corresponding values of G viz. G g 
and G^ from eqns,(2.15) and (2.17). 

4. Obtain a refined value of X from eqn.(B.4) 


where, 




f(x b ) - x b .f(x a ) 
TUT] 


(2.18) 


f (X ) = G - G_ (2.19) 

3 a 

f(X b ) = G - G b (2.20) 


The value of G in above is accurately known from 


eqn. (2.7). 

5. Check if f (X ) < e ( = 0.005). If yes stop, 

c 

else continue. 

6. Replace X Q by X^ and X^ by X^. 

7. Go back to Step 1. 
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Analysis with Harmonic Components Included 

If the harmonic components are included, then the 
expression for -"flux is 


where , 


<p = X. cos u + E X cos(u + © ) (2.21) 

1 n^l n n 

0 = Phase shift between fundamental and n 

n 

harmonic. 

The above can be written as 


m = X,cos u + Z [X cos(nu) + X sin(nu)] (2.22) 
* jl n^JL ^ ^ ns 


where 


2 2 2 

XT ' = X^ + X; o 
n nc ns 


Q ' = tan """ 1 (X/X) 

hj ns nc 


(2.23) 

(2.24) 

The Fourier series expression for f(<p) has already 
been given in eqn.(2.13). The coefficients A q and A^ are 
zero for the reasons as given for the case when harmonic 
were not considered. 

Therefore , 

f ( 9 ) = B.cos u + Z [B cos(nu) + A sin(nu)] (2.25) 
1 n£L n n 


where. 


and 


n 


B 


n 


r\ ** 

/ f(<p). cos(nu) du 


(2.26) 


% 


r\ ** 

A n = tT f f ^* sin(nu) du 


(2.27) 
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In the case with no harmonics considered , an explicit 
expression for was given in eqn.(A.S). However, it is not 
possible in this case to get explicit expressions for A n and 
B n . We must resort to numerical integration of eqns. (2.26) 
and (2.27) respectively to get these values. 


If eqns. (2.22) and (2.25) are substituted in eqn. 
(2.5) and the coefficients of various cosine and sine terms 
are equated for various values of n, then for n = 1 we will 
get the same eqns. as for the case without harmonics, i.e. 
eqns. (2.15) thru (2.17) and for n ^ 1 we will get 


where 


- n2 “ 2 X nc + nk “ X ns + B n /C th = 0 ' (2 - 28) 

~ " 2 '' )2 - nk “ X nc + V C th = 0 (2.29) 

Equations (2.28) and (2.29) are re-written as 


X„„ = 

nc 

= c 2 

A n 

+ °1 B n 

(2.30) 

X 

= c. 

A 

+ c 0 B 

(2.31) 

ns 

1 

n 

2 n 

C 1 = 


1, 

(k 2 

>0 o o , and 

+ n 2 «) 

(2.32) 

c 2 = 

■ _k_ 

no) 

41 

C 1 

(2.33) 


The iterative scheme to get X^ X n ? nd 9, is as 
follows : 

1. Assume X n (n ^ l) as zero, correspondingly and B n 
will be zero. 



28 


2. Obtain from eqn.(2.26) by numerical integration 
(trapezoidal rule has been used in the present study). 

3. Obtain a value of X^ iteratively by 'Secant* method as 
explained in analysis without harmonics. 

4. For n ^ 1, obtain A n and B n from eqns.(2.26) and (2.27) 
with above value of X^. 

5. Calculate X nc and X ng from eqns.(2.30) and (2.31) respect- 
ively. 

6. Re-compute X^ iteratively with above values of X nc , X nQ , 

A n and B n- 

7. Re-compute X and X „ from the value of X-, obtain is 

^ nc ns 1 

step 6. 

8. Repeat steps 6 and 7 until the successive values of X^, 

X and X are within a specified e (=0.005). 
n c n s 

9. Calculate X and © from eqns . (2.23) and (2.24). 

n n 

10. Calculate fundamental and harmonic components of the 
voltage by substituting, X^, X n and © n in eqfn.(2.2l) 
and then differentiating it with respect to time. 

2.3.2. NUMERICAL DATA AND RESULTS FOR STEADY STATE 

Now we give the actual data and the results for the 
field test case. As mentioned earlier, C^ includes the EMVT 
capacitance, the bus capacitance and the various bushing 
capacitances connected to the bus. The EMVT capacitance may 
be taken as 600 to 800 pf [18]. The 220 KV bus is 108 meters 
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long having A.CSR conductors of 36 mm diameter and spacing of 
4.5 m. Since there is uncertainity in the values of various 
bushing capacitances, a range of values of is taken from 
2,000 to 3,200 pf [25 j. 

The numerical values of various parameters shown in 
Fig. 2.6 are taken as 

C ± ~ 1,300 pf 

Em * J(1^220 — = 179.6 KV 

w = 314.15 rad/sec for 50 Hz system. 

The value of R has been neglected , since its effect on 
steady state solution is very small as verified by computer 
studies but not reported here, it is also reported in [34 J 

For L, the — i characteristic is as shown in Fig. 

2.8. Since the field test data is available only for the 

5 4 

slopes = 1.3x10 H and L 2 = 2.02x10 H, the slope has 
to be assumed. The results are reported for two values of 
Lg viz. 1000 H and 700 H respectively. 

The calculated results for typical values of and 
the slope Lg are given in Tables 2.1, 2.2 and 2.3. Both 
thO fundamental and the harmonic components are given in 
the Tables. It may be observed that there are three possible 
solutions. This may also be seen from Fig. 2.9, which 
depicts the condition that the source voltage e^ must be 
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FIG. 2.9 POSSIBLE SOLUTIONS FOR CIRCUIT 
SHOWN IN FIG. 2. 7. 
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TABLE - 2.1 


Possible voltage solution for fundamental and harmonic 

components 



L 3 

= 1,000 H 


S = 

2,000 pf 


Sol. 

Funda- 

3rd 

5th 

7 th 


No. 

mental 

Harmonic 

Harmonic 

Harmonic 



A 

V 3 ®3 

V 5 

% 

V 7 e 7 


1 

88.75 

0 0 

0 

0 

0 0 


2 

364.40 

51.2 2.4 

1.8 

13.2 

3.2 164.4 


3* 

263.40 

34.2 -178 

5.0 

-172 

1.6 -14.4 


* Unstable 

solution 







TABLE 

- 2.2 




Possible voltage solution 

for fundamental 

artd harmonic 


* 


components 





L 3 

= 1,000 H 


c b = 

2,500 pf 


SoI7 

funda 

3rd 


5th 

“fth 


No. 

mental Harmonic 

Harmonic 

Harmonic 



A 

V 3 0 3 

V 5 

®5 

V 7 °7 


1 

76.8 

0 0 

0 

. 0 

0 0 


2 

394.5 

50.6 2.4 

4.1 

161.0 

2.1 156 


3* 

290.9 

37.5 -172 

2.9 

-165 

2.3 -11 



* Unstable solution 
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TABLE - 2.3 

Possible voltage solution for fundamental and harmonic 

components 

L 3 = 700 H 0 b = 3,200 


Sol. 

No. 

Funda- 

mental 

3rd 

Harmonic 

V 3 9 3 

5th 

Harmonic 

V 5 ®5 

7th 

V 7 0 7 

1 

64.7 

0 

0 

0 

0 

0 

0 

2 

340.3 

62.5 

2.3 

1.5 

60.5 

3.5 

167 

3* 

273.4 

40.6 

-178 

5.1 

-171 

2.1 

-12.6 


Unstable solution. 


equal to e^ + e L (Fig. 2.7). In the first solution, the 
network appears as an inductive load.* In the second solution it 
appears as a capacitive load and the third solution is 
unstable as verified later. Regarding the high voltage (HV) 
mode, Figs. 2.10a and 2.10b show the typical voltage wave- 
forms across the EMVT for two different cases. These wave- 
forms tally well with those observed in the field (Figs. 2.2a 
and 2.2b). 

Test for stability 

Stability of the solution is tested by means of the 
conditions derived in the Appendix C (eqns. C-9 and C-10). 

The procedure for testing the stability of a particular 
solution is given below s 



VOLTAGE IN KV 


.FUNDAMENTAL 
RESULTANT 


-lOOf 


L 3 *700H 
3200 PF 

Vj (RMS)= 340.8 KV 
V 3 (RMS)s 62.5 KV 
V ( RMS) s 346.0 KV 
-THIRD HARMONIC 
i ME IN mS 

J\ 20 ms 


-30Qh 


- 4 OOf 


FIG. 2*10A sfSx PHASE VOLTAGE AT BUS ALONG 
WITH THIRD HARMONIC COMPONENT. 


FUNDAMENTAL 


RESULTANT WAVE FORM 



ttl O 
0 

S-IOO 

O 

> -20O 


La * IOOOH 


V, (RMS) * 364.4KV 
V 3 (RMS)s 51.25 KV 
V (RMS) a 368-OKV 

-THIRD HARMONIC 
if —►TIME IN ms 


-3001- 


-400r 


-500 


FIG. 2.10B \fJx PHASE VOLTAGE AT BUS ALONG 
WITH THIRD HARMONIC COMPONENT. 
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1. Calculate the angle 6 by 

6 = tan - 1 (G s /G c ) 

where G and G are calculated from eqns. (2.15) 

0 s 

and (2.16) respectively. 

2. Calculate X and Y from X , and 5 . 

00 1 0 

X = X. cos 6 

o 1 

Y q = X J _ sin 6 

3. Take the small perturbations in X q and Y q and 

obtain the partial derivatives ^ and 57 

respectively. Let they be represented by [P] 

4. Substitute [p] in eqns* (C-9) and (C-10) and 
find out whether the inequality conditions are satisfied or 
not. 

We illustrate the procedure by an example. 

X = - 168.0 is the solution for the fundamental 

component of flux which corresponds to Case I in the Table 
2.3. Since we have neglected the resistance R (Fig. 2.6), 
the angle 6 and K will be zero. 

Then eqns. (C-4) and (C-5) become 

M = -to 2 X + B 1 /C - G (2.34) 

N = -co 2 Y + Aj/C (2.35) 

Since 6 is zero, we have X q = -168.0 and Y q = Qv 
Substituting these values in eqns. (2.34) and (2.35) we get 
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M and 
0 

N . 
0 

These 

turn out 

to be 

zero since X^ and 

0 

components 

of X x 

which is 

the solution. Now, if we 

X and 
0 

Y by an 
o 

amount of 

10, we 

get 


X' 

0 

= 

X + 10 

o 


-158.0 


Y ' 
0 

= 

Y + 10 

o 

= 

10.0 

then 

dM 


M(X|, Y q ) 

- M 

0 

= -O.lOxlO 7 


dX 


X’ - X 
• o o 


and 

dN 


N(X , Y * ) 
o o 

- N 

0 

= -0.97xl0 6 


dY 


Y ’ - Y 

0 0 


Similarly 






dM 

W 

= 

0 and 

dN 

dX 

= 0 


Substituting the above in eqns. (C-9) and (C-10) 

we find that the inequalities are satisfied, and hence, 

the solution X = -168.0 is stable. 

o 

2.4. PRE-STEADY STATE (TRANSIENT) PERIOD ANALYSIS 

So far we have considered the steady state solution 
of eqn.(2.5) and have shown that two possible steady state 
stable modes are possible viz., HV mode and the LV mode. 

For studying the pre-steady state behaviour and the value 
of the initial conditions which will result in one or the 
other mode, the state space technique as given below is used. 
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Let the 

state 

variables 

be 


*1 = 

cp 



(2.36) 

y 2 = 

9 



(2.37) 

Then from eqns 

(2.5), 

(2 ,36) and 

(2.37) we get 


y l = 

y 2 



(2.38) 

y 2 = 

-k y 2 + 

f ( y x ) 

C + 

U th 

G cos(u+ 6 ) 

(2.39) 


It may be noted from eqn.(2.37) that y 2 is the 
EMVT voltage which we wish to determine. The set of first 
order differential eqns . (2.38) and (2.39) are solved by 
modified Euler method. The initial conditions y^(0) and 
y 2 (0) , which in Fig. 2.6 is the instance when the circuit 
breaker opens, are given by 


y^o) = 

- 0 ) E cos 6 
m 

(2.40) 

y 2 (°) = 

E m sin 6 
m 

(2.41) 


where s is the angle of the source voltage at the instant 

of the CB opening i.e. t = 0 

The values of the angle s between 0 to 180 degrees 

need to be considered to examine all possible cases (it 

should be obvious that the value of 5 from 180 to 360 

degrees will not give different results). The numerical 

values of the circuit parameters as given in Sec. 2.3.2 were 

used for pre-steady analysis except for the value of loss 

7 

resistance R , which is taken equal to 0.68x10 ohm 
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Tl 'd u ! »6oHOA fid u ! DllQA 


F*g. 3.11(b) Voltage across EMVT for 
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(corresponding to a loss of 2.38 KW per phase). It may be 
pointed out that this loss accounts for the EMVT , additional 
equipment at the 220 KV bus and the bus itself. It was 
observed if the losses were assumed to be somewhat larger 
(R = 0.5x10 JX.) the steady state voltage always settled down 
to the KV mode. If the losses were taken to be somewhat 
smaller (R = 1.5xlOVL) ? the steady state could not be reached 
even though the computations were carried out for a long 
time. Figures (2.11a) and (2.11b) show the waveforms of the 
voltage across the EMVT for the values of 6 equal to 40 
and 0 degree respectively. These figures may be compared 
with Figs. (2.2a) and (2.2c) obtained from the field invest- 
igations. There is a very close agreement both in magnitude 
and shape between the field and the analytical waveforms. 

The analysis was carried out for other values of 6 between 
0 to 180 degrees in step of 10 degrees. Some values of 6 
gave rise to the HV mode, while others gave rise to the LV 
mode. 

TABLE - 2.4 

Occurence 'of HV or LV mode with the angle 5 . 


S.No. 

1 

2 

3 

4 

5 

6 

7 

8 

9 

6 in degrees 

0 

10 

20 

30 

40 

50 

60 

70 

80 

Mode 

LV 

LV 

LV 

LV 

HV 

LV 

HV 

HV 

LV , 
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S.No. 

10 

11 

12 

13 

14 

15 

16 

17 

18 

6 in degrees 

90 

100 

110 

120 

130 

140 

150 

160 

170 

Mode 

HV 

LV 

LV 

HV 

HV 

LV 

> 

'3 

1 

LV 


2.5. PREVENTIVE MEASURES 

Two references [24,33] have come to our knowledge, 
where several incidents similar to that reported by us are 
given. Overvoltages due to ferro-resonance are also reported 
in [22,23,25,34,41], Some of the preventive measures 
suggested are : 

(a) To have a sequence of switching operations, such that 
the capacitive coupling between the source and the de~energi~ 
zed bus through CB grading capacitance is avoided [41]. 
However, any such recommendation which concerns with restri- 
ction on the operating procedure, may be difficult or not 
advisable to implement. 

(b) To connect a saturable reactor in parallel with the EMVT 
secondary [32]. This recommendation needs to be examined in 
terms of its technical feasibility, economics and burden on 
the EMVT. 

(c) To connect a resistance in parallel with the EMVT 
secondary [3]. As mentioned in Section 2.2 a resistance of 
9.0 JX.. in the secondary (R _ ) , which corresponds to 

O si O 

g 

0.36x10 _51L as seen from the primary side, was connected. 
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However, this did not prevent the occurence of ( the HV mode. 

In Sec. 2.4, if the value of resistance (R) to simulate the 
losses was taken as 0.68xl0uX. , then the HV mode did occur. 
While, if the value of R was taken as 0.5xloIrt_ , then the HV 
mode did not occur. Of course , these values of R are indicative 
only, since the circuit parameters like L and C etc. will also 
influence the results. The main point to note is that while 
R is of the order of 10 7 ,R^ is of the order of 10 8 

Since R and R are in parallel, unless R is reduced 
considerably, its influence on the occurence/non-occurence 
of the HV mode will not be significant. If R is reduced 

u t C 

further, the burden on the EMVT secondary may not be accept- 
able. Also the continuous heat loss due to R may be too 

O VL- 

large for a practical design, if R is installed within 

o fc- C 

the EMVT. 

* 

It is obvious that it will be better if the steady 
state overvoltages are controlled by a suitable design of 
the EMVT saturation characteristic itself rather than the 
auxiliary means listed above. To examine this situation , 
various values of the saturation slope (Fig. 2.8) were 
considered. For each value of L^, the source angle 6 was 
varied between 0 and 180 degrees to find out when the HV mode 
did occur . The result of this study is shown in Fig. 2. 12. 

It is seen from the Figure that if is greater than 775 H, 

O 

HV mode did not occur for any value of the angle 6 . This 



Slope L3 in Henrys 
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Fig, 2.12 Effect of change in L3 on the 
probability of HV mode- 
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confirms our contention that if the EMVT saturation chara- 
cteristic can be controlled, then the problem of steady state 
overvoltages along with its possible damage can be avoided. 

In order to study the problem in greater detail 
the data for values of CB grading capacitance and EMVT 
saturation characteristic for several makes was collected. 
Table 2.5 gives the values of for several makes of CBs . 

It may be noted from the Table that the range of 
is from 250 pf to 1,300 pf. The saturation characteristics 

TABLE - 2.5 

Equivalent grading capacitance ' C^' for various makes 0 f CBs. 


S.No. Make of CB 


Grading capacitance in pf .C^ 


1. Oerlikon FS9C3 


4x1,000 


= 250 


2. Magrini 245 MHM 


2x1,000 


3. BBC type DLF 


2x2,600 


1,300 


4. BHEL, India 


4x2,000 


5. Russian 


4x2,000 


of the various makes of the EMVTs are shown in Fig. 2.13, 
and as observed, they differ widely. To take an extreme 
situation it was assumed that L 2 = in Fig. 2.8. This 
means that the EMVT saturation characteristic can be approxi- 
mated by two straight line segments’ as shown in Fig. 2.14. 
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Exciting current in amps 


Fig. 2.13 EMVT saturation characteristics obtained 
from various manufacturers . 




flux 
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Current in m A 

Fig. 2.14 Typical 9-1 characteristic 
of the EMVT. 
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A number of computer studies were done to find out the cri- 
tical value of above which the HV mode did not occur for 
the values of the angle 6 from 0 to 180 degrees in steps 
of 10 degrees, keeping and fixed. Table 2.6 shows 

the results for = 1.3xl0 5 H and C ± = 1,300 pf (which 

are the values from the field test case reported earlier). 

As already mentioned, 1,300 pf is the maximum value for C^. 

A number of additional computer studies were carried out when 
the lowest value of i.e., 250 pf, was taken, keeping 

at the same value as before. 


TABLE - 2.6 

The critical value of above which the HV mode does not 
occur for the various values of break point flux 9 -^ . 


S.No. 

Break 
in p 

point flux 9 1 
.u. ± 

Critical L ^ in Henries 

1 . 

0.9 


1,200 

2 . 

1.1 


1,000 

3. 

1.3 


650 

4. 

1.5 


300 


It was observed that HV mode did not occur even if was 
reduced to a very low value of 300 H as shown in Table 2.6. 
This points to the fact that the system with higher value of 
is more prone to overvoltages. This has also been reported 
by Tsirel and Palyakov [41], 
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It is clear from Table 2.6 that if the saturation 

characteristic of the EMVT is controlled, then the over- 

voltages can be prevented. The higher the break point, lower 

can be the value of L^. However one point is clear from 

5 

the Table, compared to the slope = 1.3x10 H, the slope 

is very small (300 to 1,200 H) irrespective of the break 
point. Hence, as long as the EMVT saturation is moderate 
upto the maximum voltage expected across the EMVT, the 
problem of steady state overvoltage should not arise. Some 
manufacturers recommend [37,47] break point of 1.2 to 1.3 p.u. 
as in the event of a single line to ground fault the voltage 
on the healthy phases may rise to 1.2 to 1.3 p.u. 

2.6. CONCLUSION 

Steady state overvoltage may arise while switching a 
system which is terminated by a transformer. The actual 
damage to a 220 KV EMVT has been analysed. Subsequent to 
the damage, field investigations were performed and waveforms 
of voltages across the EMVT for the pre-steady and steady 
state periods were recorded. Analytical study using the 
principle of 'Harmonic Balance' has been employed to deter- 
mine the steady state overvoltages. It is established that 
the overvoltage occurred due to ferro-resonance phenomena. 

The validity of the method of analysis is verified by the 
fact that the steady state voltage waveforms, both in 
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magnitude and shape, are in close agreement with those 
obtained in the field. 

The pre-steady state analysis has been carried out 
using the state space technique. The validity of the analysis 
is confirmed by the fact that the waveforms obtained in the 
field and those obtained by analysis matched very well. This 
part of the analysis has also given a better insight on how 
to prevent steady state overvoltages from occurence. Based 
upon the data collected for various makes of the CBs and 
EMVTs and a number of computer studies, it has been estab- 
lished that the steady state overvoltages can be kept to a 
low value if the EMVT saturation is controlled. 
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CHAPTER III 


DYNAMIC OVERVOLTAGES 


3.1. INTRODUCTION 

In a power system, a situation can arise that due to 
sudden load rejection an unloaded EHV line is left connected 
to an over-speeding generator (s ) . This situation can lead to 
excessive overvoltages (called 'dynamic overvoltages') in the 
system. Essentially such overvoltages are the fundamental, 
frequency overvoltages with harmonic components and may 
persist for a few seconds. These voltages can become dangerous 
if the generator self-excitation takes place due to excessive 
line charging including series capacitors if any, and/or due 
to ferro-resonance which may occur because of the transformer 
saturation. These overvoltages are also influenced by the • 
synchronous generator dynamics, voltage regulator, prime-mover 
governor and network parameters. Therefore, all the above 
factors should be included for a realistic evaluation of the 
dynamic overvoltages. 

In the methods reported in literature, usually digital 
simulation is employed to predict the variation in magnitude 
and frequency of the internal voltages of the generator (s) 
and this information is used subsequently on the transient 
network analyzer (TNA) to determine dynamic overvoltages. 
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The digital simulation neglects the network transients but 
can predict self-excitation. In this Chapter the above refe- 
rred two stage computation is combined into a single digital 
simulation, which incorporates not only the network transients 
on a 3-phase basis but also the generator representation in 
detail. The prime-mover governor is also represented along 
with exciter and voltage regulator. The non-linearities due 
to the transformer and generator saturation are also included. 
It is believed that this is a more accurate approach and can 
be used to determine in a single simulation, if the phenomena 
of self-excitation and/or ferro-resonance can occur. This is 
illustrated with the help of case studies for both hydro and 
thermal systems feeding a load through a long EHV line. The 
effectiveness of the excitation system and shunt reactors 
in controlling the overvoltages is also studied. 

3.2. S YSTEM DESCRIPTION 

An isolated system consisting of a generator supplying 
load through a long EHV line has been used for the analysis 
and is shown in Fig. 3.1. When the load is rejected at the 
receiving end, the unloaded EHV line is left connected to 
the over-speeding generator. The receiving end transformer 
T 2 will be left connected to the system if it is 'low side 
switching', i.e., circuit breaker A is opened. On the other 
hand, if it is ’high side switching', i.e*, breaker B is 



5 o 



Tj « Sending end transformer 
Tj * Receiving end transformer 


Fig. 3.1 A typical EHV system 
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opened, then the transformer is not left connected to 
the system. 

The study of dynamic overvoltages upon load rejection 
will require the detailed representation of all the power 
system components shown in Fig. 3.1, which is given in the 
following sub-sections. 


3.2.1. Sy nchronous Generator 

The stator of the synchronous generator is modelled 
[28] by a dependent current source 1^ in parallel with an 
inductance as shown in Fig. 3.2. It may be noted that 

I is a 3x1 vector and L is a 3x3 matrix. 

g * g 


Such representation can handle both the symmetrical 
and unsymmetrical networks equally well. The time varying 
coupling between the rotor winding and the stator windings 
is replaced by a dependent current source and hence the 
inductance matrix L is adjusted to be constant even in the 

g 

presence of saliency. The transient or the sub-transient 
saliency is handled by introducing dummy coils on the rotor 
d or q axes. The major advantage of this circuit model of 
the synchronous generator is that the formulation of the 
system equati^teV/^ facilitated. Also, as the inductance 
matrix L ^owrT N 4^^M'.'^2 nif is a constant, this results in 
simplification of*^he * o^CTS^tions involved in the integration 
of machine equations. The ails of the mathematical formu- 
lation are given in reference [28], whereas the generator 
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winding arrangement and the final equations are given in 
Appendix - D for ready reference. 


The generator saturation is considered by reducing 
its field voltage. The reduction in field voltage E fd is 
given by 




K (E p - K ) 
e ^ J 


(3.1) 


where 

Ep = x£ Vlf (3.2) 

The coefficients and will depend upon the 

saturation characteristic of the generator. 


3.2.2. Governor system 

The block diagrams of the typical hydraulic and 
thermal governors [46] used in the study are shown in Figs. 

3.3 and 3.4 respectively. The state equations are as follows; 


For Hydraulic Governor 


*hl 

1 

= T s 

^■ u ref 

-to- (-^ 2 + 6 G h ) " G h- a] 

(3.3) 

*h2 

1 

~ T 

C- X h2 

+ 6 

(3.5) 

*h3 

1 

~ T 

w 

t Q h - 

2 <4,3 - G h>3 

(3.6) 


under the restriction that 





G h 

min 

< a 

^ ^h max 

and 

G h 

min 

< °h 

< G h max 


(3.6) 

(3.7) 
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F or T hermal Governor 


x sl 

= 

i 

T™ 

s 

8 

H 

0> 

i 

• “>/ R - G thJ 

(3.8) 

X s2 

= 

1 

' T h P 

[G th - 


(3.9) 

X s3 

= 

1 

T 

ip 

t X s2 ~ 

X s 3 j 

(3.10) 

X s4 

= 

1 

t ip 

C X s3 " 

X s4^ 

(3.11) 

under the restriction 

that 



G th 

min 

< 

X sl < 

G th max 

(3.12) 

and G th 

min 

< 

X sl < 

G th max 

(3.13) 

3.2.3. Excitation 

system 




The system considered is IEEE type - I rotating exci- 
tation system [45], whose block diagram is given in Fig. 3.5. 
The state equations for the excitation system are as follows : 


• 

E fd 

_ 1 
' T e 

[-K e .E fd + X el - S e .E fd ] 

(3.14) 

X el 

1 

“ T 

t- X el + K a (V ref * V T ' X e2> ^ 

(3.15) 

ct 



X e2 

if 

*t 

[' x e2 + K f E fd^ 

(3.16) 


under the restriction that 

V„ < X , < V u 

R min el R max 


(3.17) 




Fig *3.5 IEEE Type -I rotating excitation system. 
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3.2.4. Transfo r mer 

The equivalent circuit of the transformer is shown in 
Fig. 3.6. Its leakage inductance is represented by and 

the series resistance accounts for the copper loss. The 

magnetising inductance is represented by L which is non- 
linear due to the transformer saturation. The manufacturer 
usually provide the transformer v-i characteristic, which 
can be converted into the g -i characteristic. The cp -i 
characteristic is approximated by straight line segments, 


so that the slope L ^ is piecewise linear, as given in the 
Appendix-E. The resistances R and R account for core 

0 p 

losses. Since the core losses are voltage dependent, the 


major part of these losses is accounted for by R^. The resis- 
tance R is considered in series with L , to facilitate the 
c mt 

decay of d.c. current which exists during transient condi- 


tion. This may otherwise not decay fast enough if the source 


resistance and the transformer resistance R^ are small In 
value or are neglected. The value of R can be chosen so 
that the ratio of the magnetising reactance (X ^.) at rated 
voltage and frequency to R c is 20. Howe.ver, a lower value of 
R giving rise to a higher ratio does not change the results 
significantly. 


3.2.5. Sh un t R eact or 

The equivalent circuit of the shunt reactor is simply 
a resistance (R ) in series with inductance (L r ) . The value 
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R lt L,t 



Fig. 3 6 Equivalent circuit of the transformer. 


R, L, 


o 


— vAAflO; 

r— — O 

0 

=c 

C = 

L 


Fig. 3.7 Equivalent circuit of the transmission 

tine. 



of the L r depends upon the reactor size. The resistance 


is so chosen that the ratio of X^ r (reactor reactance at 
rated voltage and frequency) to R r is 300 to 400 A 
non-linear reactor can also be considered. In that case the 


value of will not be a constant and its representation 
will be piecewise linear similar to that of the transformer 
magnetizing inductance L 


3.2.6. Transmission Lin e 

The transmission line is represented by an equivalent 
-IS as shown in Fig. 3.7. The series resistance R^, series 
inductance and half line charging capacitance C are 

3x3 matrices for 3-phase representation. If the line is trans 
posed, then these matrices are symmetric. 

As an example, the form of L-^ is 


'si 

L ml 

Snl 

ml 

L sl 

L ml 

'ml 

L ml 

L sl 


(3.18) 


where 


L , and L , are self and mutual inductances, 
si ml 


Normally the transmission line data is given in terms 

pf its positive and zero sequence series inductances L+ and 

Lq respectively. The values of and in terms of 

and L are given by 
o 



L ml = t~ < L o - LJ (3.20) 

The form of the matrices and C is similar to 
that of the matrix L^. 

3.2.7. Load 

The load is represented by a current source whose 
magnitude and phase angle depend upon the amount of loading 
and the power factor. The load rejection is simulated by 
opening a particular phase whenever the current source value 
goes to zero in that phase. 

3.3. NETWORK STATE EQUATIONS 

Based upon the equivalent circuits given in Section 
3.2 for the various components, the system of Fig. 3.1 
is represented by the equivalent circuit shown in Fig. 3.8. 

In the Figure the R's, L's and C ! s are 3x3 matrices while 
i's, v's and cp’s are 3x1 vectors. The state variables are 
the currents i^, i ? , ig, and i ; voltages v 1 and v 2 and 
fluxes g>-j_ and <p 2 . 

The set of first order differential equations 


describing the 

system 

are 



(Lg + L 2 ) 

• 

i l 

88 V 1 

-( R 2+ R g )lj * L g X g 

(3. 

.21) 

L 3 

• 

i 2 

= V 1 

“ v 2 R 3 1 2 

(3, 

.22) 

L 5 

• 

i 3 

= V 1 

“ R 3 i 3 

(3, 

.23) 

L 6 

^4 

= V 2 

~ R^ i /, 

6 4 

(3, 

.24) 
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C 1 Vj^ = 

-ijL- x 2 “ i 3 “ f i^l^ “ P 9 lv l 

(3.25) 


C 1 *2 = 

i 2 - f 2 (<p 2 ) - i 4 + i L “ R i lv 2 

(3.26) 


9l = 

v i - R 1 f l (( Pl) 

(3.27) 


9 2 = 

v 2 - R ? f 2 ^2^ 

(3.28) 


Equations 

(3.21) to (3.28) may be written in 

concise 

form as 

[A] [X n ] 

= f ([X n ]) + [B][u(t)j 

(3.29) 

where 

[x n ] t = 

t 1 ! i 2 i 3 i 4 v l v 2 fl *2^ 

(3.30) 


[u(t)] t 


(3.31) 

faxj) 

35 

-(R 2+ R g)”V 1 

(3.32) 


- R 3 i 2 + v 1 - v 2 
" R 5 i 3 + V 1 

- ^5^4 + v 2 


0 4 4 v 

Wh - rt - Vn* 

x z J v 2 9 


- Vh - 4' f 2 (<P2>) 

V 1 ” R l^l 
V 2 - R 7 f 2 (9 2 ) 
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[A] = 

L g +L 2 

0 

0 

0 

0 

0 

0 

0 

(3.33 


0 

L 3 

0 

0 

0 

0 

0 

0 



0 

0 

5 

0 

0 

0 

0 

0 



0 

0 

0 

L 6 

0 

0 

0 

0 



0 

0 

0 

0 

C 1 

0 

0 

0 



0 

0 

0 

0 

0 

C 1 

0 
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0 

0 

0 

0 

0 

u 

0 



1 0 

0 

0 

0 

0 

0 

0 

u 


and 







1 



II 

'ca 
\ — < 

“ L g 

0 






(3 

.34) 


0 

0 









0 

0 









0 

0 









0 

0 





i 




0 

u 









0 

0 









0 

0 








It may be 

noted that 

the matrix [A] 

is block 

diagonal. 

and hence 

' , can be 

inverted eesily* Equation 

(2.29) can be 

written as 









11 

r— , 

c; 

• >< 

[A]" 

i -i 

f([X n ]) + [A] 

[B] [ 

u(t) ] 


(3.35) 

1 # „ 

or [X n J = 

?([X n j) + 

[D}[u(t) ] 




(3.36) 

where 











<p([Xj) 

, =s 

-l 

[A] 

f ([x n ]) 


| 


(3.37) 
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and [DJ = [A.]” 1 [b] (3.38) 

3 . 4 . SOLUTION PROCEDURE 

Following load rejection, the set of eqns. given 
by (3.36) is required to be solved. The components of the 
forcing function [u(t)] are the load current ij , and the 
derivative of the dependent current source 1^ as given in 
eqn. (3.31). The load current i^ is known from the load 
data and is assumed to be a constant current source. The 
expression for has been given in eqn.(D-15). It may be 
observed that it depends upon the prime-mover speed and the 
generator d and q axes currents. Also it may be noted from 
eqns. (D-6) to (D-ll) that these currents are dependent on 
the governor output P^ and the excitor voltage E fc j» The P ffi 
and can be calculated by solving the generator, governor 

and excitation system state equations. The interconnection of 
various variables for each step of integration is shown 
in Fig. 3.9. 

Before the integration procedure can be started, the 
initial values of various variables are needed. The procedure 
for the calculation of the initial conditions is given In 
Appendix-F. 
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Fig. 3.9 Interconnection of system components 
for the digital simulation. 
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3.5. CALCULATION OF DYNAMIC OVERVOLTAGES 

The calculation of dynamic overvoltages for both 
typical hydro and thermal systems are given below. The 
hydro system is found to be more prone to dynamic overvoltage 
because usually the hydro generators are smaller in size and 
are situated far away from the load centre leading to long 
EHV lines. Also the rise in speed on load rejection is more 
than for a comparable thermal system, which contributes to 
excessive overvoltages. However, the recent trend is to 
locate the thermal stations at the pit-heads and hence they 
may also become more prone to such overvoltages, especially 
when only one unit (small generation) is in operation. 

3.5.1. Case Study 1 ; Hydro Syste m 

The single line diagram of the system studied is 
shown in Fig. 3.1. The data for the various components of 
the system i.e. the generator including its excitation and 
governor system, sending and receiving end transformers along 
with their saturation characteristics, shunt reactors and the 
transmission line is given in Appendix-E. It may be observed 
from Fig, 3.1 that the load may be rejected either by opening 
of the circuit breaker A (high side switching) or circuit 
breaker B (low side switching) . The receiving end transformer 
is left connected to the system in the case of low side 
switching. Both high side and low side switching have been 
studied. It may be mentioned that If the low side switching 
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is permissible, the installation of the circuit breaker A 
could be deferred untill more than one EHV line is commi- 
ssioned. 

3. 5. 1.1. High side switching s The results for this study 
are shown in Figs. 3.10 to 3.12, where both the excitation 
system (ES) blocked and active have been considered. While 
Fig. 3.10 shows the generator terminal voltage (vy), change 
in speed (Aoj) and field voltage (Ey d ) , Figs. 3.11 and 3.12 
show the receiving end voltages on all the three phases. 

It may be noted from Fig. 3.1Ha that with ES blocked 
Vy slowly rises to 2.5 p.u. and this high voltage persists, 
indicating self-excitation. However, when ES is active, the 
maximum Vy is only about 1.45 p.u. and it settles down to 
1.0 p.u. in an oscillatory manner showing the effectiveness 
of the ES in controlling the overvoltage due to self- 
excitation. 

The rise in speed is about 40 / when ES is blocked 
and it is about 60 / when ES is active (Fig. 3.10b). This is 

v 

in general agreement with the speed rise of 40 / reported in 
the discussion of ref. [4] for a typical hydro system. The ■ 
rise in speed is less when ES is blocked as the system losses 
are higher due to persistent overvoltages. 

As shown in Fig. 3.10c, E fd is held constant at the 
initial value when ES is blocked, while it hits the 



AOO in P.U, V T in P.U. 





Fig. 3.10 Variation of Vj , AGO and Ef^ with 
time* (High side switching). 


V in P.U. V in P.U. V in P.U 


III 


w 


fstm 


ma 


0.02 0.7 0.8 1.8 

Time in seconds 

(c) Phase V voltage 


1.85 


Fig. 3.1t Receiving end voltage with ES 
active (high side switching ) . 
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Time in seconds 


(a) Phase ‘a* voltage 



Time in seconds 
(b) Phase 'b voltage 



0.02 0.7 0.8 1.8 1.85 

Time in seconds 

(c) Phase V voltage 


Fig. 3.12 Receiving end voltage with ES 
blocked (high side switching) . 
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negative ceiling and varies in ai oscillatory manner to bring 
the Vj to 1.0 p.u. when ES is active. 

Figures 3.11 and 3.12 show the receiving end voltages 
with ES active and blocked respectively. As seen from Fig. 
3.11, the maximum dynamic overvoltage (MDO) with ES active 
is 1.7 p.u., while with ES blocked it is more than 2.4 p.u. 
and still rising for the period of study (1.85 seconds). 

This further demonstrates the effectiveness of ES in cont- 
rolling the MDO. 

Effect of Shunt Reacto r 

In order to study the effect of the shunt reactor, 
it was removed from the system, with all the other data 
(including load) remaining unchanged. The receiving end 
voltages for this case are shown in Fig. 3.13. As may be 
observed the MDO is found to be 1.95 p.u. as against 
1,7 p.u. (Fig. 3.11) with shunt reactor. This demonstrates 
the effectiveness of shunt reactor in controlling the MDO. 

The shunt reactor also lowers the steady state voltage at 
the generator terminals for a given receiving end voltage 
and load at specified power factor. In the present study the 
generator terminal voltage is 1.05 p.u. and 1.12 p.u. with 
and without shunt reactor respectively. It may be mentioned 
that 1.12 p.u. voltage at the generator terminal may not be . 
acceptable for steady state system operation. 
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Time in seconds 


(a) Phase a’ voltage 



0.02 0.8 0.9 1.8 1.85 


Time in seconds 

(b) Phase V voltage 



Time in seconds 


(c) Phase V voltage 

Fig .3.13 Receiving end voltage without 
any reactor (high side switching) . 
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3. 5. 1.2. Low si de switching : To study the effect of the 
receiving end transformer left connected to the system, the 
case of ’low side switching’ was also studied. Again two 
cases are studied viz. both ES blocked and active. The 
receiving end voltages for these cases are given in Figs. 
3.14 and 3.15 respectively. 

It may be observed that the ferro-resonance takes 
place for both the cases with MDO. equal to 2.0 p.u. However, 
it is controlled when ES is active (Fig. 3. 15) and is per- 
sistent when ES is blocked (Fig. 3.14). 

It may be concluded from above that it is not 
advisable to perform the ’low side switching' as it may 
lead to ferro-resonant mode. Although the ferro-resonance 
may be eventually controlled, the MDO and the duration for 
which the overvoltage occur may cause the damage. 

3.5.2. Case Study 2 : Thermal System 

The data for a typical thermal system, where single 
line diagram is shown in Fig. 3.1, is given in Appendix - E. 
Only ’high side switching’ is considered as ’low side 
switching’ is ruled out as given in Section 3. 5. 1.2. 

For the basic system, where the EHV line length is 
280 Km, the receiving end voltages on load rejection are 
given in Fig. 3.16. The MDO in this case is 1.6 p.u. It is 
clear that the overvoltages in this case are less severe as 
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Time in seconds 
(a) Phase a* voltage 



Time in seconds 
(b) Phase b' voltage 



Time in seconds 


(c) Phase V voltage 

Fig. 3.14 Receiving end voltage with ES 
blocked (low side switching). 
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Tim® in seconds 
Ca) Phase a* voltage 



Time in seconds 
(t>) Phase b voltage 



Time in seconds 
(c) Phase V voltage 


Fig. 3.16 Receiving end voltage fora 

thermal system (high side switching). 
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compared to the hydro system of approximately "the same size 
(Fig. 3.11). However, when EHV line length is increased to 
370 Km (shunt reactor compensation kept at 60;^ of line charging) 
the MDO is increased to 2.0 p.u. for both cases (Fig. 3.17) 
and the overvoltages are persistent. 

This is because the thermal systems are more prone 
to self-excitation as they have less capability of supplying 
leading MVARS which is received through larger EHV lines. 

Also it may be observed that the ES can affect the flux along 
the generator d-axis only and hence it can not control the 
self-excitation if it occurs along the generator q-axis. 

This explains the persistent nature of overvoltage even when 
ES is active (Fig. 3.17). 

3.6. CONCLUSIONS 

A single stage digital simulation technique for 
determinate dynamic overvoltages following load rejection 
has been developed and presented in this Chapter. The various 
components of the power system have been modelled in detail 
including the network transients on a 3-phase basis. 

Both the hydro and thermal systems have been consi- 
dered. It has been shown that the excitation system is 
effective in controlling the overvoltages caused by self- 
excitation as well as ferro-resonance . The effectiveness of 
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U.7 «•' 

0 02 0,6 Time »n seconds 

(a) Phase a‘ voltage 



i nL v w : o.7 i.8 

0.02 0.6 in seconds 

(b) Phase .t> voltage 




0.7 18 

0.02 0.6 .jjjne in seconds 

(c) phase 'o' voltage 
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the shunt reactors in controlling the dynamic overvoltage is 
also demonstrated. It is also shown that the ’low side 
switching 'may lead to ferro-resonant mode and, therefore , 
should be avoided. 
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CHAPTER - IV 


TRANSIENT OVERVOLTAGES - STATE SPACE AND METHOD 
OF CHARACTERISTICS 


4.1. INTRODUCTION 

We have already studied the occurrence of the steady 
state power frequency and harmonic overvoltages in Chapter-II 
and dynamic overvoltages in Chapter-Ill. The topic which is 
yet to be covered is the study of transient overvoltages 
which can arise either due to switching or faults. An EHV 
system consists of transmission lines along with the compo- 
nents like transformers and generators etc. The transmission 
line modelling is quite complex for the calculation of 
switching transients , For this purpose several methods have 
been suggested in the past [6,7,8,9,10,11,12,20,21,29,31]. 
Perhaps the most popular method is due to Dommel [11,21,31 ] 
based upon the method of characteristics. It enables the 
transmission line to be represented by a resistive network 
along with the current sources, whose values are dependent 
upon the past history. Further, by using trapezoidal rule 
of integration, it is possible to represent the inductive and 
capacitive elements by resistances and current sources which 
depend upon the past history. Since the equivalent circuit 



for a linear network (comprising of generators, transformers 
etc.) has only L's C's, voltage and current sources, it is 
obvious that the resultant network will be resistive. The 
network admittance matrix is constant as long as the time 
step of integration is not varied. Therefore the solution of 
the network is very simple, since the admittance matrix needs 
to be triangularized only once. However, when we wish to 
include non-linear elements, then the admittance matrix is 
not constant any more, Dommel has proposed Compensation 
Method [2l] in such cases. The method works well in the case 
of one non-linear element. But if there are several non- 
linear elements, the resulting non-linear equations must be 
solved simultaneously. It is evident that the computation 
time per iteration may become excessive in such a case. For 
several non-linear elements, the state space approach may be 
better. Further, it enables the detailed generator represent- 
ation with its aroortisseur circuits, rotor dynamics and 
saturation along with voltage regulator and governor [28 J 
as already described in Chapter -III. . 

In this Chapter the two alternatives as given below 

have been, examined,, ■ 1 „ - 

(a) To interface the transmission line model based upon 

the method of characteristics with the state space model for 



82 


(b) To use state space model for the transmission line 

as well, along with rest of the system. 

4.2. SYSTEM ANALYSED 

Since the main purpose of this Chapter is to examine 
the transmission line model as interfaced with rest of the 
system, it is adequate to represent the generator by its 
simplest representation, i.e. a voltage source behind a 
constant impedance. 

The system analysed is shown in Fig. 4.1, which 
consists of a generator with unit transformer connected to an 
EHV transmission line with a shunt reactor at the receiving 
end. The representation of the generator, the unit transformer 
and the reactor is common in both alternatives 'a' and 'b' 
given above, whereas the transmission line modelling is diff- 
erent. An equivalent circuit of the system excluding the 
transmission line is shown in Fig. 4.2. As already mentioned, 
the generator is represented by a constant voltage (e g ) behind 
a constant impedance. The combined 'source* impedance 
Z (=R + jo)L ) shown in Fig. 4.2 consists of the leakage 

impedance of the transformer along with the generator 
impedance. The transformer magnetising impedance Z t is repre- 
sented by R t + jwL t and the reactor impedance Z r is represen- 
ted by R + joiL • The resistance is the circuit breaker 
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EHV Transmission 

9-HH ' l "‘ 

Generator f ranS f 0rmcr 

Fig. 4.1 The system analysed. 


a Shunt 
J reactor 



Fig. 4. 2 The equivalent circuit of the system . 
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closing resistor, which is inserted for a few milliseconds 
during the closing operations. 

The resulting differential equations for the three 
phase system are 


di 

c 

dt 

= L* 1 (-R, i c + e s - ei ) 

(4.1) 

flit 

dt " 

= L^ 1 (-R t i t + e p 

(4.2) 

di 

r 

dt 

-- L p (—ft i + g ) 
r r r m' 

(4.3) 

where e^ 

= e, + i, R, 
k km b 

(4.4) 

^m 

i c ~ x t 

(4.5) 


ir 

(4.6) 


The variable i c , i t , i r , ± km , i mk , e g , and e R 
e m are 3x1 vectors and R c , R b , R t , L c and L t are 3x3 

matrices. 


and 


4.2,1. ALTERNATIVE 'a* 

The transmission line model using the method of 
characteristics as given by Dommel [ll] is shown in Fig. 4.3. 
By using a modal transformation, the three phase* line currents 
and voltage can be converted into modal components. The 
transformation is given by 




Fig. A. 3 The transmission tine model for 
’Alternative a". 
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e phase or 1 phase = ^ e mode or Vode (4 * 7) 

where 

[T] = 1 1 1 (4.8) 

1-2 1 
11-2 

This enables the three phase transmission line to be 
treated as equivalent to three independent single phase lines. 
The resultant equations for each mode, where the super- 


script 

’prime 

’ refers to 

modal quantities, 

are given by 


i » 

x km 

= G. e£ + 

( t - t ) 

(4.9) 


i ’ 
mk 

= G * e m + 

<* - *> 

(4.10) 

where 

G * 

1/(Z + R/4) 


(4.11) 


and i.» . i\,. e’ and e* ate the modal line currents and 
km* mk* k* m 

voltages at two ends of the line* I£(t-x ) and I^(t- x ) are 
the modal current sources which depend upon the past histo- 
ries and 'x' is the modal propagation time. The Z and R 
in eqn. (4.11) are the modal surge impedance and resistance 

respectively. 

I£(t- x) and I^(t -r) are given by 

I/.(t ~c) = ^ [~ G - " i mk^ t "" T ^ ^ 

+ ( t— x ) *” lj^pj(t— x )] 


(4,12) 
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I m ( t_ ^ 

where 


— r- a - 

+ ¥ f-G-^(t-t) - (4.13) 



(4.14) 


(4.15) 


e m - 5 ^mk - yt-*)] (4.16) 

In the transmission line model presented above, there 
are four terminal quantities viz. e^, e m , i km and i mk . If any 
two of the above quantities are known, the other two can be 
calculated. In our case, first the currents i km and i k are 
obtained from eqns. (4.5) and (4.6), and then the voltages 

e. and e_ can be obtained from eqns. (4.15) and (4.16). 

km 

Now it is required to integrate eqns. (4.1) and (4.2) using 
eqns. (4.3), (4.15) and (4.16). An integration scheme using 
the Modified Euler Method! is below : 

Step 1 :■ Calculate derivatives i- c> i t and i r from eqns. (4.1) 
to (4.3). 
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Step 2 
Step 3 


Step 4 


Step 5 


Step 6 


Step 7 


Increment the time. 

Obtain the first guess for the state variables i^, 
i^. and i r with the help of derivatives calculated 
in step 1. 

Obtain the first guess for the voltages and e^ 
from eqns. (4.15) and (4.16). Note that the modal 
transformation given in eqn.(4.7) is required. 
Calculate the new values of the derivatives i c> i^ 
and i from eqns. (4.1) to (4.3). It may be noted 
that the transformation back from modal to phase 
domain is required. 

Take the average of the derivatives calculated in 
steps (l) and (5), and use this to calculate final 
value of state variables i c , i^., and i r . 

Calculate the final value of the voltages e^ and 

e' from eqns. (4.15) and (4.16). 

m 


4.2*2* ALTERNATIVE jbj 

If we also employ the state space model for the trans- 
mission line, first alternative could be to use the it-section 
model. However, it is inefficient in applications requiring 
large band width, such as simulation of electromagnetic 
transients [11,12]. The undamped natural frequencies of a 

!c 71 

loss-less, open ended line are , k = l,2 , 
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whereas, a single ^-section has an undamped natural frequency 

2 

of , which is considerably different than the line’s 

lowest natural frequency. If we partition the line into two 
te— sections, then the undamped natural frequencies are 

2f2 , 4 

^ = ~~x — and which is an improvement at the 

lowest natural frequency >but the error is large for the 

second frequency. As the line is partitioned into increas- 

I 

ingly largetf number of sections, the undamped natural fre- 
quency of t/he lowest pole approaches £ . However, the number 
of state y ariabl os is much larger and hence the model is very 

inefficient. 

The second alternative is to use the 'Fourier cosine 
series/ transformation as given by Triezenberg [44], where 
he has shown that the line model based on the representation 
of distributed voltage has correct natural frequencies with 
both ends open circuited. The model is suited to the study 
of energisation transients and any other situation where the 
line is lightly loaded. It is efficient for the calculation 
of power frequency (steady state) conditions, since it will 
require a large number of terms in the cosine series. For 
this particular situation the dual model employing the 'Fourier 
cosine series* transformation for the spatially distributed 
current is more appropriate. In view of the above, we have 
employed the model based upon the transformation of distri- 
buted voltage. 
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If v(x,t) is the three phase line voltage at a distance 
x measured from the sending end, then the 'Fourier cosine 
series ' is given by 

^[vCxjt)] = V(x,s) =“ V (s) + § Z V, (s) Cos(k7c x/d) 

° k=l K 

(4.17) 

where d is the line length and the coefficients V k (s), 
k =0,1, 2 , 3 .... are found from the integral transform 

d 

V k (s) = g- / V(x,s). Cos(k7i x/d) dx (4.18) 

o 

An approximate representation for eqn. (17) is 

V(x,s) = “V Q (s)+§ l V k (s).Cos(kTt x/d) (4.19) 

k — 1 

where *n ’ is a finite integer. 

Similar transformation can be made for the currents. 

As seen from eqn. (4.19) the voltage at any point of the line 
can b© calculated by summing the various frequency components 
[V k (s) f k« l,2,....n). Thus, the Laplace transform of the 
end point voltages e k and e m {Fig. 4.2) are given by 

X. K<t>] = E k (s> - i V 0 (s) + \ V k (s) (4.20) 

E ra (s) y ;» 0 W + | Jl ( -l )k v k ( s) 

( 4 . 21 ) 
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A transfer function which relates the terminal currents 
anc * ^mk^ s ^ ^he com P°nent voltages is derived in 
Appendix G. The structure of the line model based upon the 
eqns. (3-5), (4.20) and (4.21) is given in Fig. 4.4, where 
the value of 'n' is taken as three. The currents ^] cm ( s ) an ^ ^mk^ s ^ 
are shown in Fig. 4.4. The matrices R, L, C and G are 
3x3 matrices for a 3-phase system. 


For the sake of convenience, we define 


X 1 


i km 

■ ^mk 

(4.22) 

X 2 

= 

i km 

+ V 

(4.23) 

x 3 

as 

2 

n 

v 0 (t) 

(4.24) 

x 4 

38 

Tt 

v 2 (t) 

’ (4.25) 

x 5 

38 

2 

% 

v x (t) 

(4.26) 

x 6 

38 

SL 

% 

v 3 (t) 

(4.27) 

where x^, 

x 2’ 

X 3 , 

x 4 , , x^ , v 0 (t). 

v x (t), v 2 (t) and v 3 ( 


are all 3 x 1 vectors. 


Then from the block diagram given in Fig. 4.4 we have 

(ignoring the conductance G) 

* 3 - J 0 ' 1 n 


(4.28) 
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LC i‘ 4 + RC i 4 + (§2) 2 x 4 = | [R xj + l x x ] (4.29) 

2 

LC x 5 + RC x 5 + (f) x 5 = | [R x 2 + L x 2 J (4.30) 

and 

LC + RC x 6 + (§5)% = | [R x 2 + L x 2 J (4.31) 

Equation (4.23) to(4.3l) can be written as a set of state 


equations as follows : 
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where 


b, 

1 

as 

( 2 ) 

2 

L' 

- 1 C 1 


(4.33) 

b 2 

= 

(2) 

l" 

-1 Q-l 

R 

(4.34) 

X 

= 

*4 




(4.35) 

x 8 

= 

*5 




(4.36) 

x 9 

= 

*6 




(4.37) 

In 

the 

set 

of 

eqns . 

(4.32) 

the variables x^ and x 2 


their derivatives are known as they depend upon network 
state variables as defined by eqns. (4*5), (4.6), (4.22) 
and (4.23). Now the set of eqns. (4.32) is solved by 
'Modified Euler Method 1 (which has been already described) 
along with the eqns. (4.1) and (4.2). The line voltages e k 
and o m are given by 

e k “ x 3 + x 4 + x 5 + Xj (4. 38) 

% " X 3 + X 4 - *5 - *6 (4 ' 39) 

4*3. COMPARISON OF THE ALTERNATIVES 

Both the alternatives are compared with the help of 
two examples given as follows : 
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Example - I 

The system given in ref. [6] is considered first. 

The data for the system is given in Appendix H. The line 
energization voltage transients obtained by both the alter- 
natives are shown in Fig. 4.5. It may be noted that the 
voltage for phase 'a* is in close agreement for both the 
alternatives (and with that reported in ref. [6]). For phases 
'b' and ’ c’ also the voltage waveforms are in close agree- 
ment except that alternative 'b' give rise to superimposed 
high frequency oscillations. This may be attributed to the 
absence of the combined source impedance Z c and truncation 
of the 'Fourier Series'. 

Example - II 

The second example is taken from ref. [27], which 
includes the source impedance Z c and the transformer mag- 
netising impedance Z^ . The trapped charge and non-simultaneous 
breaker pole closing is also considered. The data for this 
example is given in Appendix H. The receiving end voltage 
waveforms obtained without and with are shown in 
Figs. 4*6 and 4.7 respectively. There is a close agreement 
between the waveforms obtained by both the alternatives 
(and with that reported in ref . [27]). The superimposed 
high frequency oscillations shown earlier in Fig. 4.5, 



Voltage in R U. Voltage in 











Voltage in P-U- Voltage in P-U. Voltage in PU. 
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Alternative a’ 



Fig. 4.6 Receiving end voltages on line energisation 
without closing resistance (Example-11) 









Voltage in PU- Voltage in RU. Voltage in P-U- 



(c) Phase V voltage 


Fig. 4.7 Receiving end voltages on line energistation 
with closing resistance (Example-!!) 
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while using the Alternative ’b' disappeared in this case. 
This is attributed to the presence of the combined source 
impedance Z . 

A comparative statement of these two alternatives is 
given in Table 4.1. The Alternative ? a* has been found to be 
superior in terms of the computation time. The computation 
time is large in Alternative 'b ! because of the larger 
computation time required per step and also the smaller step 
size. The larger (doubled) step size resulted in numerical 
instability for both the alternatives. 

Table - 4.1 

Comparative Statement for Alternative 'a* and ’b' 


S.No. Study 

Step size 

Total computation 


in msec. 

time in percent 

1. Alternative ’a* 

0.05 

100 percent 

2. Alternative * to ’ 

0.025 

272 percent 


4.4. CONCLUSION 

For determining the transient overvoltages due to 
switching, the transmission line can be modelled either 
using the method of characteristics or the state space 
approach. However, when the system consists of several non- 
linear elements, barring transmission lines, the state space 
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approach is most desirable for the rest of the network. Thus , 
the two alternatives which need to be considered are (a) to 
interface the transmission line model based upon the method 
of characteristics with the state space model for rest of 
the system, and (b) to use the state space model also for 
the transmission line along with the rest of the system. 

It is shown that it is possible to employ either of the two 
alternatives. However, the Alternative * b' has two draw- 
backs. If the source impedance is neglected, the resulting 
voltage waveforms may have high frequency oscillations 
superimposed on them. It is also shown that the total 
computation time for Alternative ’b'is approximately ,;three 
times than that for Alternative ’a’. Thus it is concluded 
that the Alternative ’a’ is superior to Alternative ’b 1 . 
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CHAPTER - V 

TRANSIENT OVERVOLTAGE DUE TO FAULTS 


5.1. INTRODUCTION 

In. Chapter IV, we have presented some modelling 
aspects for determination of the transient overvoltages in 
EHV systems. The conclusion was that if the EHV system has 
several non-linear elements which include generator and 
transformer saturation, the method of characteristics for 
the transmission line and the state space approach for the 
rest of the system should be used. In the majority of the 
EHV systems, the predominant transient overvoltages are 
those due to line energisation/4s-energisation [19,26]. Such 
overvoltages can be controlled with the help of circuit 
breaker closing resistors. When the overvoltages due to 
line energisation/ re-energisation are strongly controlled, 
the overvoltages due to fault inception and subsequent 
load rejection may become more important [15,16,17]. Fur- 
ther, if single pole switching is employed, special attention 
must be paid to the overvoltages on one of the unfaulted 
phases during a single line to ground fault. This is true, 
since a ground fault on one phase could cause a flashover 
from another phase to ground, thereby converting a single 
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line to ground fault into a double line to ground fault. 
The possibility of such an occurrences greatly diminishes 
the expected benefits of single pole switching. 

The modelling of the synchronous generator for the 
determination of the above-mentioned overvoltages should 
be detailed since there can be fundamental frequency and 
natural frequency (of the system) components besides the 
harmonic components resulting from unbalanced currents 
flowing in the armature of the synchronous generator, 
especially if it is of the salient pole type [lj. 

A number of investigators [15,16,38] have used TNA 
to study the above-mentioned overvoltages. Some of the 
inherent drawbacks of the TNA, particularly the inability 
to have detailed generator representation, have already 
been pointed out in Sec. 3.1. A mathematical model for the 
digital simulation of single line to ground fault has been 
given in ref. [17]. The method is based upon the Laplace 
Transform technique and the fault is simulated by applying 
a voltage source equal and opposite in magnitude to the 
pre-fault voltage at the fault point and then making the 
use of superposition principle. Since superposition is 
only applicable to linear systems, this approach cannot 
be used for modelling when non-linear elements are present. 
A method simulating the line to ground fault using the 
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method of characteristics, which does not require super- 
position, has been developed and presented in this Chapter. 

5.2. SYSTEM ANALYSIS 

For the purpose of analysis, the system shown in 
Fig. 5.1 is considered. Here, the power flow is assumed 
from source A towards source B. Note that this figure is 
identical to Fig. 3.1 except that the load is represented 
by a voltage source. This is done because at the time of 
fault inception on the EHV line, there will also be a 
feed from the load side. An equivalent circuit for the 
system is given in Fig. 5.2, where the source A (generator) 
is represented in details as already explained in Chapter- 
Ill. The source B (load) is represented by a constant 
voltage source along with an impedance as determined by 
the short circuit MVA. A single line to ground fault is 
assumed to occur at point ' f 1 on the line as shown in 
Fig. 5.2. The fault simulation technique is given in 
Appendix - I. 

The state equations describing the ,.system during 


fault on the EHV 

line are 

•T JJ* , 


(L g + L 2 ) i x 

= e k - (R 2 + Rg)i x - I 

• 

■'g I g 

(5.1) 

• 

L 5 ^3 

e k “ R 5 1 3 


(5.2) 

^6 "*"4 ~ 

e m " R 6 1 4 


(5.3) 
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Tj s Sanding and transformer 
Tj • Receiving end transformer 
P and Q are circuit breakers 


Fig. 5.1 System analysed for fault overvoltages. 
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(L 4 + L s > 


0 — 0 
y s m 


( R + R s ) i s 


9i 


e k - R x f X (9i) 


$2 ~ e m ~ “7 "2' 


(5.4) 

(5.5) 

(5.6) 


- r 7 f 2 (q> 2 ) 

Equations (5.1) to (5.6) can be written in concise 


form as 


where 


[A] x 


m 


f'(x m ) + v' 


t 


[V 


tv j 


- [ig i 3 


i 4 i s 9i ^ 

C(ev-L n I a ) e 


(5.7) 


(5.8) 


e ( e — e ) e k e m^ 
k _ “g~g' -k m S m X m 


f-(x m ) = 


-(Rg + ftg) 


(5.10) 


~ ^5 ^3 

- a 6 U 


-< R 4 + Rs 5 i s 

- r x f^n) 

- R 2 f 2 (9 2 ) 


[A] 


L g +L 2 


(5.11) 


L ,+L_ 
4 s 


U 


U 
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The set of eqns. given by (5.7) can be written as 
-1 -1 

x m = [A] f'(x m ) + [A] v ! (5.12) 


or 

= 

m 

f ( X ) + V 

m 

-1 

(5.13) 

where 

f (x ) 
m 

= [A] f ' (x ) 

l j v nr 

-1 

(5.14) 

and 

V 

= [A] v ! 

(5.15) 


-1 

[A] needed in the above equations can be obtained 
easily, since [A] is block diagonal. 

If source A is also represented by a constant 
voltage source than only eqn. (5.1) will be modified and 
will be similar to eqn. (5.4). 

If we wish to consider the load rejection overvoltage 
wt)ile the fault is still present, we will assume that in 
Fig. 5.1 the circuit breaker ' P ' opens* This corresponds 
to high side switching as we have already concluded in 
Chapter - III that the low side switching is to be avoided. 

We will assume that the circuit breaker has opening resistors 
as suggested in ref. [23] for the control of load rejection 
overvoltages. For the above situation the equivalent circuit 
as given in Fig. 5.2 is still applicable except that a 
resistance (opening resistor) is inserted between buses 
(2 ) and (3 ). The state equations describing the system are 
as follows. 
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(L g + L 2 J H = e k -<P 2 + Vh * L g ig < 5 ' 16 ) 

o 

L 5 i 3 = e k “ R 5 *3 (5.17) 

L 6 1 4 e m R 6 i 4 (5.18) 

(L 4 + Lg) i s - e s ” v 2 " ^4 + R s ) S (5.19) 

$1 = e k - R x ^( 91 ) (5.20) 

<p 2 = v 2 - R ? f 2 (9 2 ) (5.21) 

where v 0 = e + R, (i ~ + i.) (5.22) 

2 m b mf 4 


It may be noted that the eqns. (5,16) to (5.21) 
are same as eqns. (5.1) to (5.6) except that the term 'e ' 
in eqns. (5.4) and (5.6) is replaced by ! v 2 ' which is 
defined by eqn. (5.22) 

[6] x m = f'(x m ) + u ' (5.23) 

t 

where [x m J ? f'(x ) ancJ are a d rea dy defined by eqns. 
(5.8), (5.10) and (5.11) respectively and 

[u’] t = [<e k - L g i g ) e k e m (e s - v 2 ) e k v 2 ) (5.24) 

The set of eqns. (5.23) can be written as 

x m = f ^ x m^ + u (5.25) 

where f(x m ) as already defined by eqn. (5.14) and 
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-1 

U = [Aj u' (5.26) 

With reference to Fig. -5.2, the circuit breaker 
opening resistance will not be present between buses 
(2) and (3) when the breaker P is completely opened. 
This implies disconnection of the load along with the 
transformer from rest of the system. In this situation 
eqns. (5.19) and (5.22)are to be ignored. 

5.3. SOLUTION PROCEDURE 

The solution of eqn. (5.13) or (5. 2 £) will require 

* 

I , e , e, and e which are components of the vector 'v ! 
g s k m 

or ’u*. The I is calculated by eqn. (D-15) as already 

y 

discussed in Sec. 3.4 and the source voltage e is given. 

V 

The voltages e, and e are obtained from eqns. (I-I3) and 
k rn 

(I-l^) . It may be noted that these equations are in modal 
domain and transformation from modal to phase domain can 
be done by eqn. (4.7). The terminal currents i^^ and i m ^ 
needed in eqns. (I-I3) and (I-llf) are given by 

= - b - H - f i c ’ , i ) - e k /R 9 (5 - 27) 

and when the 

(a) circuit breaker P poles are closed 

i mf = ~ 1 s ” x 4 ~ e n/ R 8 (5.28) 
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(b) with the opening resistor inserted between poles 

trf = - *4 - f 2 (< P 2 ) - V R 8 (5 - 29) 

(c) with the poles completely open 

i mf = - i 4 (5.30) 

The solution procedure utilizing the modified Euler 
Method for the solution of the set of eqns. (5.13) or (5.2£r) 9 
(1-13) and (l-l£f) is similar to as given in Sec. 4. A,!. 

5.4. SYSTEM STUDIES 

Based on the theory presented so far, a computor 
program was developed. To gain confidence in the program 
the system studied in ref. [15] was simulated. The system 
consists of an EHV line connected to sources at both ends, 
which were simulated by constant voltage source represent- 
ation. The line to ground fault was assumed to occur at 
middle of the EHV line. In ref. [15] two values for the 
combined source inductance 'L', which consists of the source 

i 

and transformer leakage inductance (L^ + L g in Fig. 5.2), 
have been taken. They are 0 and 40 mH respectively. It 
may be mentioned that in eqn. (5.1l), the matrix [Aj has 
a term L A + L . This term can not be taken as zero, as 
otherwise [A]” 1 can not be obtained. Therefore, we chose 
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a value of 3.2 mH compared to zero in ref. [l5j. This 
corresponds to an inductive reactance of 0.001 ohm at 50 
Hertz. Compared, to the maximum overvoltage (at fault location) 
of 2.1 and 2.0 p.u. in ref. [15] for L = 0 and L = 40 mH 
respectively , we got 'values of 2.0 and 1.9 p.u. respectively. 

Further, studies were made on the system shown in 
Fig. 5.1, whose data is given in Appendix E except for the 
source B (load) , which is given in Appendix J. The gener- 
ator is assumed to be hydro because it is expected to 
give more overvoltages due to saliency [l]. Both fault 
inception and subsequent load rejection overvoltages hsive been 
considered and are presented below. 

5.4.1. FAULT INCEPTION OVERVOLTAGES 

Three studies with source A consisting of one or 
two generators in parallel, with line to ground fault 
located at the sending end or middle of the line and detailed 
(Case - I) as well as constant voltage source (Case - II) 
representation for the source A (generator) were conducted. 

The results obtained are summarized in Table 5.1, 
which shows the maximum voltage along with time of its 
occurrence at the sending end, receiving end and fault 
location. The maximum of these overvoltages (MOV) for each 
study for both Cases I and II have been underlined.' The 
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T ABLE - 5.1 

Overvoltages on healthy phases due to fault inception on 

Phase ’a’* 


Study Study details 
No. 


q Maximum voltage in p.u. 

vr (Time of occurrence in m.second) 

* Sending Receiving Fault 

end end location 

Phase Phase Phase Phase Phase Phase 
' b ’ l c ! 'b r ’c ? ’b 1 'b' 


I 


Two generators 
operating in 
parallel. Fault 
at the sending 
end of the line 


I 1.61 —2.08 1.35 -1.34 - 

(4.68) (2.52) (5.04) (1.44) 

II -1.39 - 1.90 1.23 -1.46 - 

(15.12) (2.52 (4.32) (11.88) 


II 


Two generators I 
operating In 
parallel. Fault 
at the middle of jj 
the line 


1.59 -1.94 -1.68 -1.48 1.5 -1.62 

(4.68) (1.44) (14.76) (0.72) (4.68) (.1.80) 

1.58 - 1.76 -1.39 -1.48 -1.39 1.44 

(5.4) (0.72) (14.76) (0.72) (15.12) .(11 .16) 


One generator 
operation. Fault 
jjj at the middle 
of the line. 


I -1.63 - 1.87 1.56 -1.45 -1.47 -1.53 

(14.76) (1.44) (2.88) (2.88) (6.12) (1.8) 

II -1.63 -1.72 1.45 -1.34 -1.31 -1.33 


(14.76) (0.72) (7.2) (4.32) (14.76) (1.8) 


* The fault is created when the voltage on phase 'a' at the 
fault point ift passing through its maximum. 
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high lights of these studies are as follows : 

(1) It may be observed that Case I (detailed generator 
representation) always gave rise to higher MOV compared to 
Case II (constant voltage source representation). This can 
be primarily attributed to harmonic overvoltages because of 
the unbalanced currents flowing through the armature of the 
generator, which are ignored in Case II. 

(2) It is also observed that most severe case is when the 
fault occures at the sending end. This can be attributed 
to the fact that in this case the effect of saliency is 
maximum at the generator terminals, because the external 
line reactance which reduces the percentage difference 
between the direct and quadrature axis reactances is zero. 

(3) The two generators operating in parallel i.e. high 
source MVA, give rise to higher overvoltages as compared 
to only one generator in operation. This may be explained 
as follows. The fault inception results into the initiation 
of travelling waves, which get reflected back and forth at 
the line ends. 'With high source MVA, the source impedance 
will be lower giving rise to higher reflection coefficient 
and consequent higher overvoltage. 
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For Study I, the three phase voltage waveforms at the 
sending end (where MOV Occurred) are shown in Figs. 5.3 and 
5.4 for cases I and II respectively. Similarly for Studies 
II and III, Oases I and II, Figs. 5.5 through 5.8 show the 
phase c voltage waveforms at the sending end where the 
MOV occurred. 

5.4.2. FAULT OVERVOLTAGES FOLLOWING LOAD REJECTION 

The case of load rejection overvoltages following a 
fault on the EHV line is very important [23]. In Fig. 5.1, 
it is simulated by opening of the circuit breaker P while 
the fault is still on the line. In view of the results 
obtained in Sec. 5.4.1, the detailed representation for the 
generator (Source A) has been employed as it will give higher 
overvoltage compared with the constant voltage source repre- 
sentation. The results of this study are given in Table 5.2. 

It may be observed that the MOV in this case is 2.5 p.u. 

TABLE - 5.2 

Overvoltage due to load rejection following fault 

' " Fault "To cation in ~™ Maximum overvoltages (MOV) 

Study percent of line length (The letters in parenthesis 
No. from sending end. indi cat e the phase) 

Without R b With”* (=405 ) 


I 

50 

h- 4 

• 

o> 

(b) 

1.7 

(b) 

II 

80 

2.4 

(b) 

1.9 

(b) 

III 

90 

2.5 

(b) 

2.2 

(b) 


2 . 01 - » 


Time in seconds 


Fig. 5.3 Fault inception voltage waveform 
at the sending end of the tine for 
Case- 1 and Study-! ♦ 



Fig. 5 *4 Fault inception voltage waveform 
at the sending end of the line 
for Case -II and Study-! . 





Voltage in P.U. Voltage in ku. 
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Voltage in P.U. Voltage in P.U. 
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Fig, 5. 7 Fault inception voltage waveform 
on phase *c’ at sending end of 
the line for Case- 1 and Study-Ill. 



Fig. 5. 8 Fault inception voltage waveform 
on phase c at sending end of the 
line for Case -II and Study-Ill . 
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m to load rejection on phase' b 
the line without CB opening 


—III . 



ue to toad rejection on phased 
the tine with CB opening 
ly-ISI . 
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compared with 2.08 p.u. for the fault inception case and 
1.8 p.u. for the case of load rejection without any fault 
(Chapter -III). These overvoltages were reduced to 2.2 p.u. 
with circuit breaker opening resistor, demonstrating its 
usefulness . 

These overvoltages were found to be maximum at the 
receiving end irrespective of fault location. As seen in 
Table 5.2, the maximum overvoltage occurred with the fault 
closer to the receiving end. 

For Study -III, the receiving end phase ' b' voltage 
waveforms, where MOV occurred, with and without have been 
given in Figs. 5.9 and 5.10 respectively. 

5-5. CONCLUSIONS 

The transient overvoltages due to fault inception 
and subsequent load rejection have been studied. A digital 
simulation technique, based on the method of characteristics, 
has been presented for the simulation of the line faults. 

The technique does not require the principle of super- 
position and is, therefore, applicable even if the system 
has non-linear elements. 

It is shown that the detailed generator representation 
compared with constant voltage source representation is 
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necessary as otherwise the overvoltages may be lower by 
about 10 percent. The load throw-off overvoltages following 
fault could be more severe compared with fault inception 
or load throw-off without any fault. They can be, however, 
controlled by circuit breaker opening resistors. 
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CHAPTER - VI 

CONCLUSIONS 


6.1. GENERAL 

This thesis has been aimed at developing improved or 
new methods for the calculation of steady state, dynamic and 
transient overvoltages in EHV systems. This is important 
for the proper determination of the insulation level. 

We give below a review of the significant work done 
and also scope for further work. 

6.2. REVIEW OF THE SIGNIFICANT WORK 

Chapter II deals with the determination of steady r 
state overvoltages. The failure of an electromagnetic voltage 
transformer in a practical system, when a circuit breaker 
fitted with grading capacitors was used for switching, .has 
been analysed. The analytical studies were carried out using 
the principle of 'Harmonic Balance', which gave excellent 
agreement with results obtained in the field. The above 
technique is superior compared to the describing function 
and state space techniques employed earlier to analyse the 
problem. Based upon the data collected for several makes of 
the transformers and circuit breakers, it has been shown 
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that the steady state overvoltages can be kept to a low- 
value if the transformer saturation characteristic is 
controlled. 

Chapter III deals with the determination of dynamic 
overvoltages following load rejection. A single stage digital 
simulation method has been developed for this purpose. The 
simulation includes the network representation on a three 
phase basis and detailed generator representation along with 
its excitation system and prime mover governor. The generator 
and transformer saturation is also included. It is believed 
to be more accurate and detailed than those reported in the 
literature. It is capable of predicting the generator self- 
excitation and/or ferroresonance phenomena in a single 
simulation. The effect of excitation system and prime mover 
governor parameters, transformer saturation and shunt 
reactors etc., on controlling the dynamic overvoltages can 
be studied. A number of typical studies for both hydro and 
thermal^ systems , which demonstrate the effectiveness of the 
excitation system and shunt reactors in controlling the 
overvoltages as given. It is also shown that the receiving 
end transformer should be disconnected on load rejection. 

Chapter IV deals with modelling aspects for the 
determination of transient overvoltages when several non- 
linear elements are present in the EHV system. It is shown 



122 


that in such a case it is best to model the transmission 
line based upon the method of characteristics and rest of the 
system using state space model. 

Chapter V deals with the determination of over- 
voltages due to fault inception and load throw-off subsequent 
to the fault occurrence. A technique for the simulation of 
the fault using the method of characteristics has been deve- 
loped. This is superior compared to the older methods which 
employed superposition principle for fault simulation, as 
in that case we must linearise the non-linear elements 
present in the EHV system. The fault inception overvoltages 
have been calculated both for simplified and detailed gener- 
ator representation. The work reported in the literature 
does not include detailed generator representation. It is 
shown that the simplified generator representation leads to 
optimistic results. 

6.3. SCOPE FOR FURTHER WORK 

While determining the steady state overvoltages in 

D> 

Chapter II It was found out that the occurrence of high 
voltage or low voltage mode depends on the initial conditions. 
This particular part of the analysis was carried out using 
standard state space techniques. A solution had to be obtained 
for each particular initial condition. This is obviously 
a laborious process especially when one wishes to change' some 
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of the circuit parameters. A phase plane analysis should 
be attempted. 

Thyrister controlled static shunt compensators have 
lately become popular. Their role in the control of dynamic 
overvoltages should be investigated. 

We have not considered series compensated EHV trans- 
mission lines while calculating the various overvoltages. 

This aspect should be studied further. 
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APPENDIX - A 

FO URIER EXPRESSI ON FO R f(g>) 

In Fig. A-l, it is illustrated how f (cp ) can be obt- 
ained from cp-i characteristic. Analytical value of f ( 9 ) is 
calculated piecewise. 

From 0 to u. 


f ( 9 ) 


9l <P 2 ~ $1 * cos u “ ?2 

_™ -f T ' ' + f 

^1 h 2 l 3 


X cos u 


K i + ~T. 


(A-l) 


where 


^1 ~ L L 

x L 1 l 2 


^ L 2~~ L 1 ^ + 9 ?( L 3 ~ L 2 ^ 

L 2 L 3 


(A-2) 


From to 


9 


f(<p) = t-** + 


X COS U — 


fl 


T y , X COS u 

= K 2 — c — 


(A-3) 


where 


K 


~ ) 

2 . L2 


(A-4) 


From to tc/2 


w v X cos u 

1 = f ( 9 ) = l7~~ 

1 


(A-5) 



cos u 



FIG. A.1 DERIVATION OF f (*) FROM ♦ 
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From tc/ 2 < u < it, the derivation is similar to above. 
The Fourier coefficient is given by 


1 


71 


B 1 ~ 2 % 


f f(cp) cos u du 


(A-6) 


Substituting the values of f(<p) from the eqns.(A-l) 
through (A-5) and similar equations for ti/ 2 < u < n in eqn, 
(A-6) 9 we have 

B, = § [/^(K, + ^f22_U) cos u du + f V + 

1 3 u, 2 


o 




TC— U 


n X COS u , , n / rr . A CU5 U\ 

cos u du + j — . — cos u du + J ■ -) 


1 


X cos U' 


u. 


TC-U,- 


cos u 


Tt 

du + / (-Kj, + 


X cos u 


■) cos u du ( A-7 ) 


ir-u 


By simplification 


B. 


2 rX 

n 


sin 2u 


v biu ^u, Y 1 x 

[£- (u x + — 2 — + £“* (u 2 -u x ) + ^(sin 2u 2 - sin 2^) 

3 2 


v sin 2u 0 

+ ™ (| - u 2 ) - — ^ 2- + 2 (K jl -K 2 ) sin u x + K 2 sin u 2 ] 

1 (A— 8) 


U, 


u. 


cos 


-i la 

X 


cos 


-1 £l 

X 


(A-9) 

(A-10) 


where 
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APPENDIX - B 


SECANT METHOD 

We wish to solve the non-linear equation 

f(X) = 0 (B-l) 

In Newton's method, we employ the iterative 


scheme 


X i + 1 


X 


f (x ± ) 
1 " F'TxTJ 


where 


f(X. ) - f (X. ,) 

f'(x.) = 1 - 1 --— 


X. - X. i 
1 1-1 


(B-2) 


(B-3) 


Substituting (B-3) in (B-2), we get 

X i-1 f <V " X i f(X i-l } 


X i + 1 


(B-4) 


mp-rrix-p 

The above iterative scheme is known as 'Secant' 


method 
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APPENDIX 


CONDITION OF STABILITY 


Let the solution of the equation 


5>* + km + 


9 


IM 

^th 


= G cos(cot) 


(C-l) 


for the fundamental component be 

cp = X(t) cos cot + Y(t) sin cot (C-2) 

where 

f (9 ) = cos cot + A^ sin cot (C-3) 

In the steady state, the amplitudes X(t) and Y(t) in 
eqn.(C-2) are constant. From the above equations under the 
assumption that X(t) and Y(t) vary slowly we can get [3j 

= M(x,y) = -w 2 x + kcoy + ^ - G = 0 


dY (t) 

dt 


2 

N(x,Y) = -CO Y + kcox + 


(C-4) 
0 (C-5) 


Let the steady state solution for x anc * Y X 0 

and Y . If we consider small variations a and p from the 
o 

solutions X and Y , these deviations approach zero with 
0 o 

increase of the time t, if the system is stable. 

From eqns. (C-4) and (C-5) 


da 

dt 


dM 

ax 


a + 


dM 

dY” p 


(C— 6) 
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dN 


a + 


dN ~ 

W P 


(C- 7) 


The solutions of these simultaneous equations have 
the form exp(Xt) where X is, determined by characteristic 
equation 


dM v 

dX " A 

dN 

dX 


dM 

dY 

dN v 

dY “ A 


0 


(C-8) 


The variations a and j3 will approach zero with 
increase in time, provided the real part of ,A is negative* 
This is given by Hurwitz criterian as 


dM dN 
dX “ dY 


(C~9) 


and 


dM dN dM dN 
dl 57 “ dT dl 


> 0 


(C-10) 
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APPENDIX - D 


GENERATOR EQUIVALENT CIRCUIT 


The generator winding arrangement is shown in Fig. 
D-l. The dependent current source model of the generator is 
shown in Fig. 3.2. 

The dependent current source vector I having 3-phase 
components is expressed as 

Ig = ^ £ + I q S (D-l) 

The transformation matrices £ and £ are defined as 

c T = 1/2/3 [cos 9 r cos(9 r - 2 tt/ 3) cos(9 r + 2u/3)J (D-2) 

s T = V"2/3[ sin 0^ sin(9 r ~ 2%/3) sin(9 r + 2^/3)] (D-3) 


The generator terminal current i^ can also be split in its 
d-q components as 


i 


d 



(D— 4) 


i 


q 



(D-5) 


Now the final equations describing the machine can be written. 
Here <p refers to flux linkage and the subscripts refer to 


winding shown in Fig. D-l. 

Z d = C 1 + C 2 9h 

*q = C 3 + C 4 'Pk 

9f ” a i ( Pf + a 2 cp h + b l E fd + b 2 1 d 


(D-6) 
(D— 7 ) 
(D-8) 




132 


9 h = a 3 9 f + a 4 9 h + b 3 i d 

Vg = a 5 fg + a 6 fk + b 5 i q 

?k = a 7 Tg + a 3 fk + b 6 1 q 

The electro-mechanical power is 


(D-9) 

(D-10) 

(D-ll) 


-x V (i , I - i I.) 
d d q q d 


(D-12) 


where x^' is sub-transient reactance. Further, 


(D-13) 


m <- D6 + p m ~ P e> 


(D-14) 


The various coefficients shown in above equations are 
dependent on self and mutual reactances of the windings shown 
in Fig. D-l. Details are available in ref. [28 ], 

cp^ , cp h , 9g> » and 6 are the state variables. 

The derivative of I , which is required for the 
calculation of forcing function, is obtained by differentiat- 
ing eqn.(D-l) with respect to time. 


(I, + 0)1 ) C + (I - W I.)s.‘ 


(D-15) 
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APPENDIX - E 

SYSTEM DATA FOR DYNAMIC OVERVOLTAGE STUDY 


Data for both the cases of study, viz., hydro and 
thermal systems are given below. All the data is on 100 
MVA base. 

a) Case Study I - Hydro System 



Generator parameters 





Rating 173.7 MVA, 

11 KV, 

50 Hz 

x d 


0.607 j 

X i = 

0.135 

x d 

= 

0.0853 ; 

X q = 

0.336 

X q 

— 

o 

— » 

x" = 

q 

0.1325 

X 

o 

=s 

0.065 ; 

R 

a 

0.00158 

T i 

ao 

= 

8.84 ; 

T If — 

do 

0.03 

T * 
qo 

= 

O 

9 

T i? — 

qo 

0.03 

H 

c 

=s 

4.56 ; 

D 

c 

0.0 

N 

= 

0.035 i 

*2 " 

7.358 


K 3 = 0.8 

Line parameters 

One three phase transposed single circuit, 400 KV , 
260 Km long line. The line parameters are as follows : 
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z 1 = 0.0048 + j 0.0517 p.u. 

z q = 0.0375 + j 0.1351 p.u. 

y x = 3 1.475 

y Q = j 0.985 p.u. 

These symmetrical component parameters are converted 
into self and mutual impedances in the digital programme. 


Tra nsformer parameters 

A typical non-linear saturation characteristic is 
taken from ref. [23 J and is shown in Fig. E-l. This chara- 
cteristic is voltage versus current which is converted into 
flux versus current. The characteristic is taken as piecewise 
linear by four straight line segments as shown in Fig. E-rl. 


Sending end transformer 
Rating t 180 MV A , 


R, 


R 


0.0 p.u. 
277.7 p.u. 
4.0 p.u. 


11 KV/400 KV 

; = 0.0833 p.u, 


(accounts for a loss of 0.65 MW ) 
R, = 4.0 p.u. (accounts for a loss, of 0.09 MW) 

Magnetising current at no load = 1.25 percent 

With reference to Fig. E-l where magnetising chara- 
cteristic is shown 

5 9b 

S S, = 


9 a 

9 C 


0.0032 p.u. 
0.00382 p.u. 
0.0075 


0.00352 p.u. 

0.2525 

0.00245 


0.0011 
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Receiving 

end Transformer 


Rating 

; 450 MV A , 

400 KV/ 200 KV 

H 

0.00 p.u. 

; X 4 = 0.0334 p.u 


Rg = 800.0 p.u. (Accounts for a loss of 0.562 MV) 

R^ = 0.79 p.u. (Accounts for a loss of 0.350 MV) 

The nature of magnetising characteristic is same as that of 
sending end transformer. The break point fluxes and the 
slopes are 


fa 

= 

0.0032 p.u. 

; 9 b 

sr 

0.00353 

f 0 

rr 

0.00382 p.u. 

O 

9 

A 

= 

0.101 

S 2 

= 

0.003 


S 3 

= 

0.00098 

S 4 

= 

0.00044 





Voltage 

regulator parameters 




K a 

5= 

, 400 

o 

9 

T a 

= 

0.05 sec 

K e 

= 

-0.17 

• 

t 

T e 


0.95 sec 

K f 

= 

0.04 

o 

t 

T f 

= 

1.0 sec 

S 

= 

0.00027 

9 

S b 

= 

1.304 

V D 

R max 

= ' 

3.5 

» 

9 

V R 

min 

= - 3.5 

Hydraulic 

Governor parameters 




T r 

= 

5.0 sec 

© 

9 

T 

w 


= 1.0 sec 

T s 

= 

0.2 sec 

o 

9 

a 


= 0.05 

G, * 

h mm 

= 

0.0 

© 

9 

°h 

max 

= 1.0 

a 

s= 

0.3 

o 

9 

G h 

min 

= . -5.0 

^h max 

= 

5.0 
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Shunt Reactor at Receiving end 
Rating : 65 MVAR 


X 

Load 


1.34 p.u. 


R = 0.0077 p.u. 


1.5 p.u. at J V | = 1.0 and power fact-dr = 0.95 log, 

b) Case Study II - Thermal System 
Generator parameters (Data on 100 MVA base) 

Rating : 250 MVA. ; 50 Hz 


x d 

= 

0.8463 

p.u. 

© 

9 

x d 

= 

0.1178 p.u. 

x d' 

ss 

0.0816 

p.u. 

© 

9 

X q 

SS 

0.8468 p.u. 


ss 

0.1178 

p#u.« 

o 

9 

X * 

= 

0.0816 p.u. 

X 1 

rr 

0.0652 

p.u. 

o 

9 

x o 

ss 

0.044 p.u. 

T ' 
do 

= 

7.0 sec 

9 

T * t 

do 

ss 

0.1 sec 

T' 

3 - 

0.5 sec 

1 

T ?? 

= 

0.05 sec 

qo 





qo 



R 

=r 

0.00072 

p.u. 

© 

f 

H 

= 

2.2 

B. 




c 



D 

c 

= 

0.0 


© 

9 

K i 

= 

0.035 

K 2 

ss 

7.358 


o 

9 

K 3 

s= 

0.8 


One three phase, single circuit, transposed, 400 KV 

and 280 Km long line. The parameters are as follows % 

z ± = 0.0052 + j 0.0575 p.u. 

z = 0.0452 + j 0.1760 p.u. 

o 

= j 1.495 p.u. 
y = j 1.03 p.u. 
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Sending En d Transformer parameters 
Rating : 250 MV A, 11 KV/400 KV 


R, 


R. 


R- 


0.0 


X r 


= 0.0639 


= 722.00 p.u. (Corresponding to a loss of 
= 2.88 p.u. (Corresponding to a loss of 


Details of magnetising characteristics (Fig. E-l ) 


<Pc 

S 2 

S 4 


= 0.0032 

= 0.00382 
= 0.0054 

= 0.000792 



Receiving end Transformer parameters 


0.00351 

0.1818 

0.001764 


Not considered. 


Voltage Regulator parameters 


T e 

= 


0.7 

sec 

e 

f 

\ - 

1.0 


T 

a 

= 


0.5 

sec 

0 

9 

K 

a 

0.12 


h 

= 


2.0 

sec 

© 

* 

h = 

200 


S a 

= 


0.0 


0 

9 

s b - 

0.0 


R max 


4.0 


* V 
’ V R 

min ~ 

-4.0 


Thermal 

qovernor 

parameters 




T sl 


= 

0.425 sec 

9 

R 

0 

= 

0.04 

T hp 


= 

0.2 

sec 

o 

9 

F hp 

= 

0.62 

T. 

ip 


= 

4.0 

sec 

0 

9 

F ip 

zz 

0.38 

T 

IP 


= 

100.0 

sec 

© 

9 

> 

= 

0.0 


min 

= 

-1.18 

0 

9 

^th max 

= 

0.59 

G th 

min 

ss 

0.0 


© 

9 

^th max 


0.85 


.137 MW) 

.082 m) 
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Shunt Reactors at Bot h ends 

Rating : 45 MVAR calculated at voltage 1.05 p.u. 

X = 2.45 ; ' R = 0.008 

Load 

2.0 p.u. at j V ! = 0.95 and power factor = 0.95 lag. 
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APPENDIX - F 

CALCULATION OF INITIAL CONDITIONS 


Network Initial conditions 

Since the load conditions are known/ the various 
voltages and currents can be calculated by performing the 
load flow. The entire network is symmetrical and balanced. 
The voltages and currents are calculated in phasor form and 
then the instantaneous values for each phase are cal- 

culated. The load current is chosen as a reference phasor. 


Generator Initial conditions 

The vector diagram for the generator is shown in 
Fig. F-l. The terminal voltage v^ and current i^ is known 
in phasor- form from the initial conditions of the network 
variables. The field voltage and the angle 5 * are 

calculated with the help of the following 
C 


vt 


V T - i l- (R a +JX q ) 


(F-l) 

6 = J_ C vt - n/2 (F-2) 

E fd = l C vt - 1 d ^ x d - x q^ (F-3) 

Now the initial value of other variables are 
given by ; 


I _ =E 

x f 


•fd 


df 


(F— 4) 
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9f 


x f . I f + x df .i d 

(F-5) 

Vh 

= 

x dh * X d + x f h * I f 

(F-6) 


=r 

X q' g * i q 

(F-7) 

9k 

= 

X q ° K * i q 

(F— 8) 

• 

6 


0.0 

(F-9) 

Governor Initial conditions 



The 

equations for calculating initial 

conditions 

are given below for both hydro and thermal governors. 

Hydro Governor 



*hi 

= G h 

(F-rlO) 


X h2 

= 6 . G h 

(F-ll) 


X h3 

ll 

fOjOJ 

• 

(F-12) 

where 

G h 

100 p 

P base * m 

(F-13) 

Under steady 

state conditions P = P^ 

m e 


Thermal 

Governor 



X sl 

- Q th 

(F-14) 


X s2 

= a th 

(F— 15) 


X s3 

= G th 

(F-16) 

t 


X s4 

- G t h 

* 

(F-17) 

where 

Excitor 

r 100 m 

th P, * F, + F. -f F. 

base hp ip lp 

Initial conditions 

(F-18) 


The 

equations are as follows s 



X el 

= 0.0 

(F-19) 


X a o 

= (S Q + K ) E- rt 

e e id 

(F-20) 
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APPENDIX - G 

DERIVATION OF THE TRANSFER FUNCTION RELATING TERMINAL 
CURRENTS TO THE COMPONENT VOLTAGES OF 'FOURIER COSINE SERIES' 


The basic equations for a transmission line in 
Laplace transform domain are : 

Z(s) I(x,s) = - V(x,s) (G-l) 

Y (s) V(x,s) = - (G-2) 

Differentiating the eqn. (O-l) and substituting in eqn. 

(A-2) we get 

,2 

Z(s) Y(s) V(x,s) = V(x,s) (G-3) 

dx^ 

If both sides of eqn. (G-3) are transformed by eqn. (4.16) 
then we get 

Z(s) Y(s) V k (s) = -(kix/d)V k (s) ^ E k (s) 

+ f^> k ^ E m^) (G “ 4) 

Finally, substituting eqn. (G-l) for the derivatives in 
eqn. (G-4) . 

V k (s) = ~ [Z(s) Y(s) + (k7c/d) 2 Uj _1 Z(s) [l km (s)-(-D k I mk (s)] 

(G-5) 
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APPENDIX - H 


DATA FOR SYS TEM CONSIDERE D I N CHAPTER IV 


Example -__I 

AXl the data as given in ref. [6] is on a 345 KV and 

100 MVA base. 
r = 0.00918 p.u. 


r = 0.0509 p.u. 

o 


X = 0.099 p.u. 


Y = 0.96 p.u. 

o 


iC x = 0.258 p.u. 

✓ = 1.64 p.u. 

* 4 . 

, ^^7 anri the transformer magnetising 

The source impedance Z c and t 

impedance Z t are neglected. 

Circuit breaker closing resistor R fe and the trapped 

charges were taken as zero. All the circuit breaker poles ^ 

closed simultaneous. The source voltage e g which is 

dered to be balanced three phase, is given below. 


e = 

1.0 sin(377t) 

s 

1.0 sin(377t - 271/3) 


1.0 sin(377t + 2%/ 3) 

HxamEiG-t 

II 

The data as given in ref. 


and 100 MVA base. 



145 


The line data 


R 1 


0.0048 

p.u. 

R 

0 

0.0375 

p.u. 

X 1 

=s 

0.0517 

p.u. 

x o = 

0.1551 

p.u* 

*1 


1.4706 

p.u. 

' Y o = 

0.9804 

p.u* 

data 






R c 

= 

0.000 

p.u. 

X c = 

0.1313 

p.u. 

R t 

ss 

0.000 

p.u. 

x t = 

80.0 

p.u. 

R r 

as 

0.0077 

p.u . 

X r = 

1.54 

p.u. 

R. 

5 z 

0.25 

p.u. 





ole closing span in msec. - 0 , 2, 4 

_ _1 o i.o, 1.0 

rapped charge in p.u. 

1.148 Sin(314t + 9 ) 1 

1.148 Sin (314t - 2 ti /3 + 9 ) 

1.148 Sin(314t + 2n;/3 + 9> ' 

here 9 = 280 degrees (this corresponds to the 

ngle which gave rise to maximum overvoltage). 
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APPENDIX - I 

S IMULATION OF T HE TRAN SMI SSION LINE FAULTS 


With reference to Fig. JL2, the fault is assumed to 
occur at point ' f' on the transmission line connected 
between terminals k and m. The transmission line is section- 
alised into sections ! a' and 'b* at the fault point as shown 
in the figure. Each section is represented as given in 
ref. [2lJ. The equations in modal domain (primed variables) 
for both the sections are as follows : 


where 


1 kf 

= G a * 

0 J 
k 

4 * 

IjJ (t - x) 

(I-D 

x fk 

= G a * 

e f 

+ 

I^ a (t -*c) 

(1-2) 

i fm 

= V 

e f 

4 - 

I| b (t - t) 

(1-3) 

i mf 

= G b . 

e ii 

4 - 


(1-3) 



+~Rj4 


(1-5) 




( 1 - 6 ) 


1+h. 


1-h 


[-G .e l(t-x)- ilu(t-x)] + - 5 -^ [- G a- e k (t " x ) 


I^t-Tc) = — 2 L-ui a .e^v t-T>- T ~2 




(1-7) 
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I.T (t-*) = -*-H.-G a - e k 


q-G = .e'(t-x) -ifcf (t- 1) 3 ^[- S a e f (t ' ,) 




( 1 - 8 ) 




' 1+h 




" i fm (t - T)] 


L i-G b . e f( t -'c) 


(1-9) 


I’(t-x) = 

m 


1+lV 

— ^[-G, e' 
2 b m 






- x: f (t-T) ] 


( 1 - 10 ) 


G a - R a /4 
r C+ R a /4 

Gb * V 4 


( 1 - 11 ) 


(1-12) 


By re-arranging the eqns. (I-D and (1-4) we get 


i- [it* - x i (t - T) 3 

a 


i- Ti’ r - I' (t - x) ] 
TT^ L x mf m' 


(1-13) 


(1-14) 


G b L mf m 


Equations (1-13) and (1-14) are similar to eqns. 
(4.15) and (4.16) and can be solved as discussed in Chapter 
IV. The determination of the terms I£(t-0 and I m (t- t 
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will require the past histories of voltage and currents 
ijcf* i fk’ if m and “Snf* The curre hts i^ f and i^ f are known 
by eqns. (5.27) to (5.30) whereas the currents i^ and i£ m 
are calculated by eqns. (1-2) and (1-3). The expression for 
will be different under different system conditions and 
is derived below. 


(a) No fault condition s 

It may be observed from Fig. 5.2, that under no fault 
condition we have 


’f . 


ifk 


" ifm 


(1-15) 

(1-16) 


From eqns. (I-l), (If 2) , (1-15) and (1-16) we have 




G + G. 
a b 


(1-17) 


(b) Line to ground fault 

A. line to ground fault is created by forcing the 
voltage on phase 'a’ at the fault point to zero. The faults 
^ phases ’b* and ' c’ will be obviously zero. 


e 


fa 


i 


fb 



(1-18) • 
(1-19) 


and 


0 
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From eqns . (4.7), the modal quantities 1,2,3 


given by 


e fl 


e f2 

1! 

CO|H* 

e f3 

L 


111 
1 -1 o 

1 0 -1 


e f a 
e fb 
e fc 


From eqns. (1-18) and (1-20) we get 


e fa = e fl + e f 2 + e f3 ° 


The fault current i f in terms of its modal 


nents is 


i fi 

I i f 2 

i f3 


1 1 1 


^fa 

1 

H* 

l 

H 

O 


i fb 

3 

1 0 -1 


he 

-r — 


From eqns. (1-19) and (1-22) we get 

fa y 


ifl = 42 " i f3 lfa/3 


At the fault point 'f' 

i h + 4m + 4 - ° 

From eqns. (1-2), (X-3) and (1-24) we get 
-[11 + I| a (t- ') * ^b (t ' 4]/(V G b> 


are 


( 1 - 20 ) 


( 1 - 21 ) 

compo- 


( 1 - 22 ) 


(1-23) 


(1-24) 


(1-25) 
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and from eqns 


X 2- 

. (I-lt), (M3) and (p-25) 


'f 1 


- ^ [X fa (t- *) + ^b^* ~ d 4 ^ Ifa 2 (t T) 

u 4 1 x ■ * 


+ Iru (t- *0] -07 «-*fa 3 


u fb. 


^ [I fa (t-T) + I fb (t^)] C 1 " 26 ) 

4 fa * 3 


where 


D 1 = 


Dr 


D, 


(G a 2 + G b2^ ^ G a3 + 

(G a3 + ‘W^al + ‘W 

(G al + G bl )(G a2 + G b2 ) 


tnd 


D 4 = D 1 + °2 + °3 


(1-27) 

(1-28) 

(1-29) 

(1-30) 


The suffix 1,2,3 indicate the modal component of 
;he quantities. 

(c) Double li ne to qround_fault 

A double line to ground fault is simulated by 

t t-v i -nr\r\ * r * at the fcfcu.lt point 
forcing the voltage on phases b and 

i- i i will obviously 
to zero. The fault current on phase a i fa 

be zero. (Mi) 


'fb 


* = 0 
f c 


"fa 


and 


0 


(1-32) 
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Substituting the eqns. (1-31) and (1-32) in eqns. 
(1-20) and (1-22) respectively we get 

(1-33) 


5 fl 


e f 2 S f3 


e f a^ 3 


and 


if a = A fl + i f2 + i f3 

f ' 


= 0 ' 


(1-34) 


'fa 


From eqns. (1-25). (1-33) and (1-34) we get 

= - b- ^ + Ifa 2 (t " T) + Ifa3<t " T) + Ifbl(t " T) 

5 1 

(1-35) 
(1-36) 


♦ + If »3 (t - T)] 


where D, = « al + G a2 + <\ 3 + G b! + G b2 + G b3 
Re-arranging eqn. (1-25), we get 

i • = - [If a (i'" T ) + I fb Ct - T)] " (G a + Gb)-6 f 


(1-37) 
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APPENDIX -> J 

DAT A fob SOURCE B CONS lDEREQji_S^^ 

mm t m mm***® * ***** I , * > 


All data is on 


100 MVA base. 


(a) r„, 1t Inception overvoltages 

Study -JL 

Positive sequence impedance 
Zero sequence impedance 
Voltage at bus (4) 

Power transfer 


= 0.0008 + 3 0.0426 p.u. 
= 0.0015 + 3 0.065 p.u. 
= 1.0 p.u. 

- 3.0 + j I*° P.u* 


St.udv - .II 

Same as Study - I 


Study - III- 

Positive sequence impedance 
Zero sequence impedance 
Voltage at bus (4) 

Power transfer 


= 0.0016 + j 0.085 p.u. 
= 0.0030 + j 0*13 P* u * 
= 1.0 p.u. 

_ t . 5 + j 0.5 p.u. 


r„ llnw -i nn load throw-off i 

(b) Faul^over v olta ges f ol i o * 

d IIX the source a 

For all the Studies - i» l 

^ * c +nd v - I of the fault, inception, 

data is same as that o 


overvoltages. 

The GB opening 


resistance 


is taken as 400 ohms 


and is inserted for 


8 m. seconds. 
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